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Abstract 

I review the results of hadron spectroscopy calculations from lat- 
tice QCD for an intended audience of low energy hadronic physicists. 
I briefly introduce the ideas of numerical lattice QCD. The various 
systematic errors, such as the lattice spacing and volume dependence, 
in lattice QCD calculations are discussed. In addition to the discus- 
sion of the properties of ground state hadrons, I also review the small 
amount of work done on the spectroscopy of excited hadrons and the 
effect of electromagnetic fields on hadron masses. I also discuss the 
attempts to understand the physical mechanisms behind hadron mass 
splittings. 

1 INTRODUCTION 

QCD at low energies is hard to solve, perhaps too hard for mere mortals 
to solve, even when assisted with the latest supercomputers. QCD is the 
theory that describes the interactions of quarks and gluons. QCD has been 
well tested in high energy scattering experiments where perturbation theory 
is valid. However, QCD should also describe nuclear physics and the mass 
spectrum of hadrons. Hadron masses depend on the coupling {g) like M ~ 
g-1/9 hence perturbation theory can't be used to compute the masses of 
hadrons such as the proton. 

The only technique that offers any prospect of computing masses and 
matrix elements non-perturbatively, from first principles, is lattice QCD. In 
lattice QCD, QCD is transcribed to a lattice and the resulting equations are 
solved numerically on a computer. The computation of the hadron spectrum 
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using lattice QCD started in the early 80's P3I21- The modern era in lattice 
QCD calculations of the hadron spectrum started with the results of the 
GFll group jHl llj. The GFll group were the first to try to quantify the 
systematic errors in taking the continuum and infinite volume limits. 

The goal of a "numerical solution" to QCD is not some kind of weird and 
misguided reductionist quest. Our inability to solve QCD has many profound 
consequences. A major goal of particle physics is to look for evidence for 
physics beyond the standard model of particle physics. One way of doing 
this is to extract the basic parameters of the standard model and look for 
relations between them that suggest deeper structure. To test the quark 
sector of the standard model requires that matrix elements are computed 
from QCD [S]. The problem of solving QCD is symbolically summarised 
by the errors on the quark masses. For example, the allowed range on the 
strange quark mass in the particle data table |6j is 80 to 155 MeV; a range 
of almost 100%. The value of top quark mass, quoted in the particle data 
table, is 174.3 ± 5.1 GeV As the mass of the quark increases its relative 
error decreases. The dynamics of QCD becomes simpler as the mass of the 
quarks gets heavier. Wittig has reviewed the latest results for the light quark 
masses from lattice QCD [7 Irrespective of applications of solutions to QCD 
to searches for physics beyond the standard model, QCD is a fascinating 
theory in its own right. QCD does allow us to test our meagre tools for 
extracting non-perturbative physics from a field theory. 

In this review I will focus on the results from lattice gauge theory for the 
masses of the light mesons and baryons. I will not discuss flavour singlet 
mesons as these have been reviewed by Michael jHUH]. There has been much 
work on the spectroscopy of hadrons that include heavy quarks jTUl HH [12] j 
however I will not discuss this work. The treatment of heavy quarks (charm 
and bottom) on the lattice has a different set of problems and opportunities 
over those for light quarks. Although the spectroscopy of hadrons with heavy 
quarks in them can naturally be reviewed separately from light quark spec- 
troscopy, the physics of heavy hadrons does depend on the light quarks in the 
sea. In particular the hyperfine splittings are known to have an important 
dependence on the sea quarks [12]. 

Until recently, the computation of the light hadron spectrum used to 
be just a test of the calculational tools of lattice QCD. The light hadron 
spectrum was only really good for providing the quark masses and estimates 
of the systematic errors. However, the experimental program at places such 
as the Jefferson lab ^3 E] has asked for a new set of quantities from 
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lattice QCD. In particular the computation of the spectrum of the A^*'s is 
now a goal of lattice QCD calculations. 

As the aim of the review is to focus more on the results of lattice calcu- 
lations, I shall mostly treat lattice calculations as a black box that produces 
physical numbers. However, "errors are the kings" of lattice QCD calcu- 
lations because the quality and usefulness of a result usually depends on 
the size of its error bar, hence I will discuss the systematic errors in lattice 
calculations. Most of systematic errors in lattice QCD calculations can be 
understood using standard field theory techniques. I have also included an 
appendix ^ on some of the "technical tricks" that are important for lattice 
QCD insiders, but of limited interest to consumers of lattice results. How- 
ever, it is useful to know some of the jargon and issues, as they do effect the 
quality of the final results. 

There are a number of text books on lattice QCD. For example the books 
by Montvay and Munster [12], Rothe [T7], Smit [TH] and Creutz jTH] provide 
important background information. The large review articles by Gupta PU] , 
Davies ^T] and Kronfeld [21] also contain pertinent information. The annual 
lattice conference is a snap-shot of what is happening in the lattice field 
every year. The contents of the proceedings of the lattice conference have 
been put on the hep-lat archive for the past couple of years [221 123 1^ - 
The reviews of the baryon spectroscopy from lattice QCD by Bali [23] and 
Edwards [20] describe a different perspective on the field to mine. There used 
to be a plenary reviews specifically on hadron spectroscopy at the lattice 
conference |2Zl 123 123 EDI EI] ■ The subject of hadron spectroscopy has now 
been split into a number of smaller topics, such as quark masses. 

If the reader wants to play a bit with some lattice QCD code, then the 
papers by Di Pierro [^211231211, contain some exercises and pointers to source 
code. The MILC collaboration also make their code publicly available (try 
putting "MILC collaboration" into a search engine). 

2 BASIC LATTICE GAUGE THEORY 

In this section, I briefly describe the main elements of numerical lattice QCD 
calculations. Quantum Chromodynamics (QCD) is the quantum field theory 
that describes the interactions of elementary particles called quarks and glu- 
ons. The key aspect of a quantum field theory is the creation and destruction 
of particles. This type of dynamics is crucial to QCD and one of the reasons 
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that it is a hard theory to solve. 

In principle, because we know the Lagrangian for QCD, the quantum field 
theory formalism should allow us to compute any quantity. The best start- 
ing point for solving QCD on the computer is the path integral formalism. 
The problem of computing bound state properties from QCD is reduced to 
evaluating equation [TJ 



where Sp and Sq are the actions for the fermion and gauge fields respectively. 
The path integral is defined in Euclidean space for the convergence of the 
measure. 

The fields in equation[T]fiuctuate on all distance scales. The short distance 
fiuctuations need to be regulated. For computations of non-perturbative 
quantities a lattice is introduced with a lattice spacing that regulates short 
distance fiuctuations. The lattice regulator is useful both for numerical cal- 
culations, as well for formal work [32] (theorem proving), because it provides 
a specific representation of the path integral [TJ 

A four dimensional grid of space-time points is introduced. A typical 
size in lattice QCD calculations is 24^ 48. The introduction of a hyper-cubic 
lattice breaks Lorentz invariance, however this is restored as the continuum 
limit is taken. The lattice actions do have a well defined hyper-cubic sym- 
metry group [2EIEZI- The continuum QCD Lagrangian is transcribed to the 
lattice using "clever" finite difference techniques. 

In the standard lattice QCD formulation, the decision has been made to 
keep gauge invariance explicit. The quark fields are put on the sites of the 
lattice. The gauge fields connect adjacent lattice points. The connection 
between the gauge fields in lattice QCD and the fields used in perturbative 
calculations is made via: 



The gauge invariant objects are either products of gauge hnks between quark 
and anti-quark fields, or products of gauge links that form closed paths. All 
gauge invariant operators in numerical lattice QCD calculations are built 
out of such objects. For example the lattice version of the gauge action is 
constructed from simple products of links called plaquettes. 




(1) 



(2) 




(3) 
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The Wilson gauge action 

Sg = -/?(E ^Tr(f/p + Ul) - 1) (5) 

is written in terms of plaquettes. Nc is the number of colours. The action 
in El can be expanded in the lattice spacing using equation El to get the 
continuum gauge action. 

Sg = E T^<F,^Fn + 0{a') (6) 



The coupling is related to P via 



The coupling in equation [7| is known as the bare coupling. More physical 
definitions of the coupling |3H] are typically used in perturbative calculations. 
The fermion action is generically written as 

Sf = '^Mi) (8) 

where M is called the fermion operator, a lattice approximation to the Dirac 
operator. 

One approximation to the fermion operator on the lattice is the Wilson 
operator. There are many new lattice fermion actions, however the basic 
ideas can still be seen from the Wilson fermion operator. 

Sj = E('^E{V^x(7m - l)U,{x)ij.+,-^,^^{^, + l)[/t(a:)^J (9) 

X fl 

The Wilson action can be expanded in the lattice spacing to obtain the 
continuum Dirac action with lattice spacing corrections, to the Dirac La- 
grangian. The k parameter is called the hopping parameter. It is a simple 
rescaling factor that is related to the quark mass via 

«^ = -T7^ 7 (10) 

2(4 + m) ^ ^ 

at tree level. An expansion in k is not useful for light quarks, because of 
problems with convergence, however for a few specialised applications, it is 
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convenient to expand in terms of k jSH]. The fermion operator in equation |H1 
contains a term called the Wilson term that is required to remove fermion 
doubling fT^ El QB] ■ The Wilson term explicitly breaks chiral symmetry so 
equation gets renormalised. Currently, there is a lot of research effort in 
designing lattice QCD actions for fermions with better theoretical properties. 
1 briefly describe some of these developments in the appendix 1X1 

Most lattice QCD calculations obtain hadron masses from the time sliced 
correlator c(t). 

c(t) = (Ee^''O(x,t)O(0,0)t) (11) 

X 

where the average is defined in equation [TJ 

Any gauge invariant combination of quark fields and gauge links can be 
used as interpolating operators (0(x, t)) in equation [TT] An example for an 
interpolating operator in equation ^2 for the p meson would be 

t\ = iJl{x, t))^^^2{x, t)) (12) 

The interpolating operator in equation El has the same J^^ quantum num- 
bers as the p. 

The operator in equation El is local, as the quark and anti-quarks are at 
the same location. It has been found to better to use operators that build in 
some kind of "wave function" between the quark and anti-quarks. In lattice- 
Q CD-speak we talk about "smearing" the operator. An extended operator 
such as 

0(x, 0) = E /(r)V^i(x, 0))7,V2(x + r, 0)) (13) 

r 

might have a better overlap to the p meson than the local operator in equa- 
tionEl The /(r) function is a wave-function like function. The function / is 
designed to give a better signal, but the final results should be independent 
of the choice of /. Typical choices for / might be hydrogenic wave functions. 
Unfortunately equation El is not gauge invariant, hence vanishes by Elitzur's 
theorem jlU] . One way to use an operator such as equation El is to fix the 
gauge. A popular choice for spectroscopy calculations is Coulomb's gauge. 
On the lattice Coulomb's gauge is implemented by maximising 

F=Y.j:m^)+ui{x)) (14) 

X i=x 
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The gauge fixing conditions that are typically used in lattice QCD spec- 
troscopy calculations are either: Coulomb or Landau gauge. Different gauges 
are used in other types of lattice calculations. 

Non-local operators can be measured in lattice QCD calculations, 

OfeO) =^^i(a;,0))7.V'2fe + r,0)) (15) 

r 

however these do not have any use in phenomenology |4H |321 UHl HH 

because experiments are usually not sensitive to single ijjijj Fock states. 

There are also gauge invariant non-local operators that are used in cal- 
culations. 

0{x,t) =i^i{x,t)liF{x,x + r)'ilJ2{2l + L,t) (16) 

where F{x, x + r) is a product of gauge links between the quark and anti- 
quark. 

There are many different paths between the quark and anti-quark. It has 
been found useful to fuzz [IB], the gauge links by adding bended paths: 

4 

PSU{3) [cUold + Uu-bend) (17) 
1 

where Psu(3) projects onto SU(3). There are other techniques for building 
up gauge fields between the quark and anti-quark fields, based on computing 
a scalar propagator jl?! |3H| • 

The path integral in equation is evaluated using algorithms that are 
generalisations of the Monte Carlo methods used to compute low dimensional 
integrals [IH]- The physical picture for equation^Jis that a hadron is created 
at time 0, from where it propagates to the time t, where it is destroyed. This 
is shown in figure H Equation ^2 can be thought of as a meson propagator. 
Duncan et al. [20] have used the meson propagator representation to extract 
couplings in a chiral Lagrangian. 

The algorithms, usually based on importance sampling, produce sam- 
ples of the gauge fields on the lattice. Each gauge field is a snapshot of the 
vacuum. The QCD vacuum is a complicated structure. There is a community 
of people who are trying to describe the QCD vacuum in terms of objects 
such as a liquid of instantons (for example pT). The lattice QCD community 
are starting to create publicly available gauge configurations HH] • This 
is particularly important due to the high computational cost of unquenched 
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Figure 1: Physical picture of a meson propagator on the lattice. The thick 
lines are the quark propagators. 

calculations. An archive of gauge configurations has been running at NERSC 
for many years. 

The correlator c(t) is a function of the light quark propagators (M(f/(i))~^) 
averaged over the samples of the gauge fields. 

1 ^ 

<t) ~ 17 E fiMm)r\ Mm))-') (is) 

i 

The quark propagator is the inverse of the fermion operator. The 

quark propagator depends on the gauge configurations from the gauge field 
dependence of the fermion action (see El for example). 

The correlator in equation^! is essentially computed in two stages. First 
the samples of the gauge configurations are generated. These are the build- 
ing blocks of most lattice QCD calculations. The correlators for specific 
processes are computed by calculating the hadron correlator for each gauge 
configuration and then averaging over each gauge configuration. 

To better explain the idea behind computing quark propagators in lattice 
QCD calculations, it is helpful to consider the computation of the propaga- 
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tor in perturbation theory. The starting point of perturbative calculations 
is when the quarks do not interact with gluons. This corresponds to quarks 
moving in gauge potentials that are gauge transforms from the unit configu- 
ration 

A^{x) = ^ U^{x) = 1 (19) 

Under these conditions the quark propagator can be computed analytically 
from the fermion operator using Fourier transforms (see [17, for derivation). 

M"™* = —. ^- — (20) 

Eto7Msin(p^) + M(p) 

2 

M{p) = m+ -^sin^{ap^/l) (21) 

Although in principle the propagator in equation 1201 could be used as a basis 
of a perturbative expansion of equation UHl the physical masses depend on 
the coupling like M ~ e~^^^ (see section so hadron masses can not be 
computed using this approach. The quark propagator on its own is not gauge 
invariant. Quark propagators are typically built into gauge invariant hadron 
operators. However, in a fixed gauge quark propagators can be computed 
and studied. This is useful for the attempts to calculate hadron spectroscopy 
using Dyson-Schwinger equations [5^ l55j. 

In lattice QCD calculations the gauge fields have complicated space-time 
dependence so the quark propagator is inverted numerically using variants of 
the conjugate gradient algorithms jSEI- Weak coupling perturbation theory 
on the lattice is important for determining weak matrix elements, quark 
masses and the strong coupling. 

The sum over x in the time slice correlator (equation fTT|) projects onto 
a specific momentum at the sink. Traditionally for computational reasons, 
the spatial origin had to be fixed either at a point, or with a specific wave 
function distribution between the quarks. Physically, it would be clearly 
better to project out onto a specific momentum at the origin. The number of 
spatial positions at the origin is related to the cost of the calculation. There 
are new lattice techniques called "all-to-all" that can be used to compute a 
quark propagator from any point to another without spending prohibitive 
amounts of computer time [HZl 110] • 

There have been some studies of the point to point correlator jOOl IHH E2] 

(O(x,t)O(0,0)t) (22) 



9 



It was pointed out by Shuryak [021 that the correlator in equation |221 niay 
be more easily compared to experiment. Also, the correlator |221 could give 
some information on the density of states. In practise the point to point 
correlators were quite noisy. 

The physics from the time sliced correlator is extracted using a fit model JH] '■ 

c{t) = aoexp{—mot) + aiexp{—mit) + ■ ■ ■ (23) 

where mo (mi) is the ground (first excited) state mass and the dots represent 
higher excitations. There are simple corrections to equation 1221 for the finite 
size of the lattice in the time direction. In practise as recently emphasised 
by Lepage et al. jMj, fitting masses from equation |2S1 is nontrivial. The 
fit region in time has to be selected as well as the number of exponentials. 
The situation is roughly analogous to the choice of "cuts" in experimental 
particle physics. The correlators are correlated in time, hence a correlated 
should be minimised. With limited statistics it can be hard to estimate the 
underlying covariance matrix making the /dof test nontrivial [0^1 EHl- As 
ensemble sizes have increased the problems resulting from poorly estimated 
covariance matrices have decreased. 

It is usually better to measure a matrix of correlators so that variational 
techniques can be used to extract the masses and amplitudes (671 EHl UHl • 

c{t)A B = (E e^'''0(2I, t)AO(0, 0)1,) (24) 

x_ 

The correlator in equation |^ is analysed using the fit model in equation |2S1 

c{t)AB= E Xanc-'^'^Xna (25) 

N=l..No 

The matrix X is independent of time and in general has no special struc- 
ture. The masses can be extracted from equation |2H1 either by fitting or by 
reorganising the problem as a generalised eigenvalue problem [69j. Obtaining 
the matrix structure on the right hand of equation EHl is a non-trivial test 
of the multi-exponential fit. One piece of "folk wisdom" is that the largest 
mass extracted from the fit is some average over the truncated states, and 
hence is unphysical. This may explain why the excited state masses ob- 
tained from |4S] in calculation with only two basis states were higher than 
any physical state 

The main (minor) disadvantage of variational techniques is that they 
require more computer time, because additional quark propagators have to 
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be computed, if the basis functions are "smearing functions". If different 
local interpolating operators are used as basis states, as is done for studies of 
the Roper resonance (see section [7^, then there is no additional cost. The 
amount of computer time depends linearly on the number of basis states. 
Also the basis states have to be "significantly" different to gain any benefit. 
Apart from a few specialised applications ,70| the efficacy of a basis state is 
not obvious until the calculations is done. 

Although in principle excited state masses can be extracted from a multi- 
ple exponential fit, in practise this is a numerically non-trivial task, because 
of the noise in the data from the calculation. There is a physical argu- 
ment [711 [72] that explains the signal to noise ratio. The variance of the 
correlator in equation ^2 is: 

4 = ^(((O(t)O(0))^)-c(t)2 (26) 

The square of the operator will couple to two (three) pion for a meson 
(baryon). Hence the noise to signal ratio for mesons is 



c{t) 

and noise to signal ratio for baryons is 



C(t) 



The signal to noise ratio will get worse as the mass of the hadron increases. 

In figure [2] I show some data for a proton correlator. The traditional way 
to plot a correlator is to use the effective mass plot. 

When only the first term in equation [221 dominates, then equation [23 is 
a constant. There is a simple generalisation of equation [221 for periodic 
boundary conditions in time. 

In the past the strategy used to be fit far enough out in time so that 
only one exponential contributed to equation [221 One disadvantage of doing 
this is that the noise increases at larger times, hence the errors on the final 
physical parameters become larger jMj . 



11 



Many of the most interesting questions in hadronic physics involve the 
hadrons that are not the ground state with a given set of quantum numbers, 
hence novel techniques that can extract masses of excited states are very 
important. 

Recently, there has been a new set of tools developed to extract physical 
parameters from the correlator c{t) in equation 1221 The new techniques are 
all based on the spectral representation of the correlator in equation ITTl 

c{t) = / p{s)e-''ds (30) 
Jo 

The fit model in equation 1231 corresponds to a spectral density of 

p{s) = ao5(s — mo) + ai6{s — mi) (31) 

Spectral densities are rarely a sum of poles, except when the particles can't 
decay such as in the Nc ^ oo limit. Shifman reviews some of the physics in 
the spectral density [74 . As the energy increases it is more realistic to repre- 
sent the spectral function as a continuous function. Leinweber investigated a 
sum rule inspired approach to studying hadron correlators [72] on the lattice. 
Allton and Capitani [7^] generalised the sum rule analysis to mesons. The 
sum rule inspired spectral density is 

Pis) = ^S{s - M) + e{s - .o)p^^^(^) (32) 
The operator product expansion is used to obtain the spectral density. 

no n-1 

pOPE^^^ = E u—TVyCnO^m,, (: qq :)) (33) 

n=l l" 

where rriq is the quark mass and : gg : is the quark condensate. The C„ 
factors in equation 1221 are obtained by putting the OPE expression for the 
quark propagator in El The final functions that are fit to the lattice data 
involve exponentials multiplied by polynomials in 1/t (where t is the time). 
These kind of methods can also help to understand some of the systematic 
errors in sum rule calculations [7^. This hybrid lattice sum rule analysis is 
not in wide use in the lattice community. See [77] for recent developments. 
There has been an interesting calculation of the mass of the charm quark 
using sum rule ideas and lattice QCD data [75] . 
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Figure 2: Effective mass plot for proton correlators. Unquenched data from 
UKQCD [73 with k=0.1350 , /3 = 5.2. 
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There is currently a lot of research into using the maximum entropy 
method to extract masses from lattice QCD calculations. Here I follow the 
description of the method used by CP-PACS [7^] . There are other variants of 
the maximum entropy method in use. The maximum entropy approach (HU] 
is based on Bayes theorem. 

P{F I DH) oc P{D I FH)P{F \ H) (34) 

In equation EH F is the spectral function, D is the data, and H is the prior 
information such as /(w) > 0, and P{F \ DH) is the conditional probability 
of getting F given the data D and priors H. The replacement of the fit model 
in equation 031 is 

Cf{t) = / K{w,t)f{w)dw (35) 



where the kernel K is defined by 

ir(w,t) = e-'"* + e-"'(^-*)) (36) 
with T the number of time slices in the time direction. 



P{D I FH) = —e-^ (37) 
L = ^EW^O-c/(tO)Aj'(c(t,)-c/(t,)) (38) 

where Zl is a known normalisation constant and is the number of time 
slices to use. The correlation matrix is defined by 

1 ^ 

A, = _ E(c(t.) - c"(t.))(c(t,) - c"(t,)) (39) 

1 ^ 

cW = T7Ec"(t.) (40) 

n=l 

where N is the number of configurations and c"(t) is the correlator for the 
n-th configuration. Equation EH is the standard expression for the likelihood 
used in the least square analysis. 

The prior probability depends on a function m and a parameter a. 

P{F I Hma) = tt^c"^ (41) 
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with Zgia) a known constant. The entropy S is defined by: 



E[/.-m.-/.log(A)] (42) 



i—i ' 

The a parameter, that determines the weight between the entropy and L, is 
obtained by a statistical procedure. The physics is built into the m function. 
CP-PACS build some constraints from perturbation theory into the function 
m. The required function / is obtained from the condition 

^i^Sjf) - L) _ 

JJ \f=U- U (4dj 

In principle the maximum entropy method offers the prospect of comput- 
ing the masses of excited states in a more stable way than fitting multi- 
exponentials to data. It will take some time to gain more experience with 
these types of methods. I review some of the results from these techniques 
in sections and [THJ 

A more pragmatic approach to constraining the fit parameters was pro- 
posed by Lepage et al. [HI] . The basic idea is to constrain the parameters of 
the fit model into physically motivated ranges. The standard is modified 
to the xlug- 

Xlug = X' + X'""' (44) 
with the prior xlrior defined by 



\2 



Xaug = 22 + 22 (45) 

n '-'an n '-'mn 

The priors constrain the fit parameters to m„ = ± am„ and an = dn ± , 
The derivation of equation Incomes from Bayes theorem plus the assumption 
that the distribution of the prior distribution is Gaussian The idea is to 
"teach" the fitting algorithm what the reasonable ranges for the parameters 
are. Lepage et al. jM] do test the sensitivity of the final results against the 
priors. 

It is difficult to understand why certain things are done in lattice calcu- 
lations without an appreciation of the computational costs of lattice calcula- 
tions. The SESAM collaboration [HI] estimated that the number of fioating 
point operations (Nfiop) needed for nj =2 full QCD calculations as: 

Nfiop cx ' (46) 

a arupi 
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for a box size of L, lattice spacing a, and Ngampie is the number of sample 
of the gauge fields in equation [TTJ A flop is a floating point operation such 
as multiplication. With appropriate normalisation factors equation EUl shows 
how "big" a computer is required and how long the programs should run 
for. Equation shows that it is easy to reduce the statistical errors on the 
calculation as they go like ~ , ^ , but more expensive to change the 

■\/ -^sample 

lattice spacing or quark mass. 

In some sense equation HHl (or some variant of it) is the most important 
equation in numerical lattice QCD. To half the value of the pion mass used 
in the calculations requires essentially a computer that is seven times faster. 
Equation 1321 is not a hard physical limit. Improved algorithms or techniques 
may be cheaper. There are disagreements between the various collaborations 
that do dynamical calculations as to the exact cost of the simulations. In 
particular, a formulation of quarks on the lattice called improved staggered 
fermions |H2] has a much less pessimistic computational cost. The reason for 
this is not currently understood. 

3 SYSTEMATIC ERRORS 

The big selling point of lattice QCD calculations is that the systematic errors 
can in principle be controlled. To appreciate lattice QCD calculations and to 
correctly use the results, the inherent systematic errors must be understood. 

For computational reasons, an individual lattice calculation is done at 
a finite lattice spacing in a box of finite size with quarks that are heavy 
relative to the physical quarks. The final results from several calculations are 
then extrapolated to the continuum and infinite volume limit. It is clearly 
important that the functional forms of the extrapolations are understood. 
Most of the systematic errors in lattice QCD calculations can be understood 
in terms of effective field theories [23]. A particularly good discussion of 
the systematic errors from lattice QCD calculations in light hadron mass 
calculations is the review by Gottlieb |27.. 

The MILC collaboration make an amusing observation pertinent to error 
analysis ^83j. In the heavy quark limit, the ratio of the nucleon mass to 
rho mass is 3/2, requiring one flop to calculate. This is within 30 % of the 
physical number (1.22). This suggests for a meaningful comparison at the 
3a level, the error bars should be at the 10% level. 
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(a) (b) (c) 

Figure 3: Different contributions to the pion correlator from Chen jS3] 

3.1 UNQUENCHING 

The high computational cost of the fermion determinant led to development 
of quenched QCD, where the dynamics of the determinant is not included in 
equation^ hence the dynamics of the sea quarks is omitted. Until recently 
the majority of lattice QCD calculations were done in quenched QCD. I will 
call a lattice calculation unquenched, when the dynamics of the sea quarks 
are included. 

The integration over the quark fields in equation can be done exactly 
using Grassmann integration. The determinant is nonlocal and the cause 
of most of the computational expense in equation 0^1 The determinant de- 
scribes the dynamics of the sea quarks. Quenched QCD can be thought of 
as corresponding to using infinitely heavy sea quarks. 

Chen |H1] has made some interesting observations on the connection be- 
tween quenched QCD and the large Nc limit of QCD. Figure El shows some 
graphs of the pion two point function. In unquenched QCD all three dia- 
grams in figure 01 contribute to the two point function. In quenched QCD, 
only diagrams (a) and (c) contribute. In the large (number of colours) 
limit only graphs of type (a) in figure (21 contribute. This argument sug- 
gests that in the large Nc limit quenched QCD and unquenched QCD should 
agree, hence for the real world SU{3) case, quenched QCD and unquenched 
QCD should differ by 0{1/Nc) ~ 30%. Chen IH firms up this heuteristic 
argument by power counting factors of Nc and discusses the effect of chiral 
logs. This analysis suggests that the quenching error should be roughly 30%, 
unless there is some cancellation such that the leading 0{1/Nc) corrections 
cancel. 

Quenched QCD is not a consistent theory, because omitting the fermion 
loops causes problems with unitarity. Bardeen jH2] et al. have shown that 
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there is a problem with the non-singlet O"*"^ correlator in quenched QCD. 
The problem can be understood using quenched chiral perturbation theory. 
The non-singlet 0+"'" propagator contains an intermediate state of t]' — tt. The 
removal of fermion loops in quenched QCD has a big effect on the 1]' propaga- 
tor. The result is that a ghost state contributes to the scalar correlator, that 
makes the expression in equation |2S1 inappropriate to extract masses from 
the calculation. Bardeen et al. IHSI predict that the ghost state will make the 
Oo mass increase as the quark mass is reduced below a certain point. This 
behaviour was observed by Weingarten and Lee [HHI for small box sizes (L 
< 1.6 fm). The negative scalar correlator was also seen by DeGrand [HI] in 
a study of point to point correlators. Damgaard et al. jHZj also discuss the 
scalar correlator in quenched QCD. 

One major problem with quenched QCD is that it does not suppress zero 
eigenvalues of the fermion operator |88J. A quark propagator is the inverse of 
the fermion operator, so eigenvalues of the fermion operator that are zero, or 
close to zero, cause problems with the calculation of the quark propagator. 
In unquenched QCD, gauge configurations that produce zero eigenvalues in 
the fermion operator are suppressed by the determinant in the measure. 
Gauge configurations in quenched QCD that produce an eigenvalue spectrum 
that cause problems for the computation of the propagator are known as 
"exceptional configurations". Zero modes of the fermion operator can be 
caused by topology structures in the gauge configuration. The problem with 
exceptional configurations get worse as the quark mass is reduced. The new 
class of actions, described in the appendix El that have better chiral symmetry 
properties, do not have problems with exceptional configurations 

Please note that this section should not be taken as an apology for 
quenched QCD. As computers and algorithms get faster the parameters of 
unquenched lattice QCD are getting "closer" to their physical values. Hope- 
fully quenched QCD calculations will fade away. 

3.2 LATTICE SPACING ERRORS 

Lattice QCD calculations produce results in units of the lattice spacing. One 
experimental number must be used to calculate the lattice spacing from: 

a = amfatt/^expt (47) 

As the lattice spacing goes to zero any choice of m^^j should produce the 
same lattice spacing - this is known scaling. Unfortunately, no calculations 
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are in this regime yet. The recent unquenched calculations by the MILC 
collaboration [HHl EOl EI] may be close. 

Popular choices to set the scale are the mass of the rho, mass splitting 
between the S and P wave mesons in charmonium, and a quantity defined 
from the lattice potential called tq. The quantity tq [121 HH] is defined in 
terms of the static potential V{r) measured on the lattice. 

'^'"^ ^-^^ ^^^^ 

Many potential |l92j models predict tq ~ 0.5 fm. The value of tq can not 
be measured experimentally, but is "easy" to measure on the lattice. The 
value of tq is a modern generalisation of the string tension. Although it may 
seem a httle strange to use ro to calculate the lattice spacing, when it is not 
directly known from experiment. There are problems with all methods to set 
the lattice spacing. For example, to set the scale from the mass of the rho 
meson requires a long extrapolation in light quark mass. Also it is not clear 
how to deal with the decay width of the rho meson in Euclidean space. 

The physics results from lattice calculations should be independent of the 
lattice spacing. A new lattice spacing is obtained by running at a different 
value of the coupling in equation [7| In principle the dependence of quantities 
on the coupling can be determined from renormalisation group equations: 

P{go) = «^^?o(«) (49) 

The renormalisation group equations can be solved to give the dependence 
of the lattice spacing on the coupling [HI] . 

a = i(^7o^7o)-^^/(^^«^e(-^/(2-o^o))(l + 0(^,0^)) (50) 

The e~^/^^'^°^0'* term in equation EU] prevents any weak coupling expansion 
converging for masses. 

Equation EU] is not often used in lattice QCD calculations. The bare 
coupling qq does not produce very convergent series. If quantities, such as 
Wilson loops are computed in perturbation theory and from numerical lattice 
calculations the agreement between the two methods is very poor. Typically 
couplings defined in terms of more "physical" quantities, such as the plaque- 
tte are used in lattice perturbative calculations 
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Figure 4: The lattice spacing determined from the mass of the rho meson 
as a function of the couphng (/5 = p-) ^] 

AUton [HH] tried to use variants of EOlwith some of the improved couphngs. 
He also tried to model the effect of the irrelevant operators. An example 
of Allton's result for the lattice spacing determined from the p mass as a 
function of the inverse coupling is in figure 01 

There are corrections to equation EDI from lattice spacing errors. For ex- 
ample the Wilson fermion action in equation El differs from the continuum 
Dirac action by errors that are 0{a). The lattice spacing corrections can be 
written in terms of operators in a Lagrangian. These operators are known 
as irrelevant. Ratios of dimensional quantities are extrapolated to the con- 
tinuum limit using a simple polynomial. 

^ = ^ + - + 0(a^) (51) 
The improvement program discussed in appendix ^ designs fermion actions 
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to reduce the lattice spacing dependence of ratios of dimensional quantities. 
When computationally feasible, calculations are done at (least) three lattice 
spacings and the results are extrapolated to the continuum. This was the 
strategy of the large scale GFll ^ and CP-PACS calculations jHE]- 

It is very important to know the exact functional form of the lattice 
spacing dependence of the lattice results for a reliable extrapolation to the 
continuum. The quantum field theory nature of the field renormalises the 
correction polynomial. There are potentially 0{a"'g"^) errors. This is par- 
ticularly a problem for states with a mass that is comparable to the inverse 
lattice spacing. There is physical argument that the Compton wavelength of 
a hadron should be greater than the lattice spacing. This is a problem for 
heavy quarks such as charm and bottom, that is usually solved by the use of 
effective field theories [TUl 1111 1^ . The masses of the excited light hadrons 
are large relative to the inverse lattice spacing, so there may be problems 
with the continuum extrapolations. 

There have been a number of cases where problems with continuum ex- 
trapolations have been found. For example Morningstar and Peardon jHT] 
found that the mass of the 0"'""'" glueball had a very strange dependence on 
the lattice spacing. Morningstar and Peardon had to modify the gauge 
action to obtain results that allowed a controlled continuum extrapolation. 
The ALPHA collaboration jHlI discuss the problems of extrapolating the 
renormalisation constant associated with the operator corresponding to the 
average momentum of non-singlet parton densities to the continuum limit, 
when the exact lattice spacing dependence is not known. 

In principle the formalism of lattice gauge theory does not put a restric- 
tion on the size of the lattice spacing used. The computational costs of lat- 
tice calculations are much lower at larger lattice spacings (see equation I46p. 
However, there may be a minimum lattice spacing, set by the length scale 
of the important physics, under which the lattice calculations become unre- 
hable [Tlin] . 

There has been a lot of work to validate the instanton liquid model on the 
lattice |lUlj . Instantons are semi-classical objects in the gauge fields. The 
instanton liquid model models the gauge dynamics with a collection of in- 
stantons of different sizes. Some lattice studies claim to have determined that 
there is a peak in the distribution of the size of the instantons between 0.2 to 
0.3 fm in quenched QCD |1021 11031 1104j . If the above estimates are correct, 
then lattice spacings of at least 0.2 fm would be required to correctly include 
the dynamics of the instanton liquid on the lattice. However, determining 
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the instanton content of a gauge configuration is non-trivial, so estimates of 
size distributions are controversial. At least one group claims to see evidence 
in gauge configurations against the instanton liquid model [ 1U51 llU6j . 

There have been a few calculations of the hadron spectrum on lattices 
with lattice spacings as coarse as 0.4 fm P71 11U7| HQS. 10^ IllOj . There were 
no problems reported in these coarse lattice calculations, however only a few 
quantities were calculated. 

Another complication for unquenched calculations is that the lattice spac- 
ing depends on the value of the sea quark mass . as well as the coupling. 
This is shown in figure El The dependence of the lattice spacing on the sea 
quark mass is a complication, because for example the physical box size now 
depends on the quark mass. 

Some groups prefer to tune the input parameters in their calculations so 
that the lattice spacing [5^lll2j is fixed with varying sea quark mass. Other 
groups prefer to work with a fixed bare coupling |29|j. 

3.3 QUARK MASS DEPENDENCE 

The cost of lattice QCD calculations (see equation HU)) forces the calculations 
to be done at unphysical large quark masses. The results of lattice calcula- 
tions are extrapolated to physical quark masses using some functional form 
for the quark mass dependence. The dependence on the light quark masses 
is motivated from effective field theories such as chiral perturbation theory. 
The modern view is that the lattice QCD calculations do not need to be done 
at exactly the masses of the physical up and down quarks, but the results 
can be matched onto an effective theory |113j . In this section I briefly discuss 
the chiral extrapolation fit models used to extrapolate hadron masses from 
lattice QCD data to the physical points. 

The relationship between effective theories and lattice QCD is symbiotic, 
because the results from lattice QCD calculations can also test effective field 
theory methods. The extrapolation of lattice QCD data with light quark 
masses is currently a "hot" and controversial topic in the lattice QCD com- 
munity. The controversy is over the different ways to improve the conver- 
gence of the effective theory. There was a discussion panel on the different 
perspectives on using effective field theories to analyse lattice QCD data, at 
the lattice 2002 conference |114j . In section |31 I show that the masses of the 
quarks used in lattice calculations are getting lighter, that will help alleviate 
many problems. 
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Figure 5: Lattice spacing dependence on the sea quark mass l/^sea from 



In the past chiral perturbation theory provided functional forms that 
could used to analyse the quark mass dependence of lattice data. In prac- 
tise, the extrapolations of lattice data were mostly done assuming a linear 
dependence on the quark mass |3j , The mass dependence from effective field 
theory was essentially treated like a black box. However, a number of high 
profile studies |115j have shown that a deeper understanding of the physics 
behind the expressions from chiral perturbation theory is required. The sub- 
ject of effective Lagrangians for hadronic physics is very large. The reviews 
by Georgi [116J, Lepage [117^ and Kaplan |118j contain introductions to the 
physical ideas behind effective field theory calculations. Scherer |119j reviews 
chiral perturbation theory for mesons and baryons. 

The chiral perturbation theory for the lowest pseudoscalar particles is the 
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perhaps the most "well defined" theory. The effective field theories for the 
vectors and baryon particles are perhaps more subtle with concerns about 
convergence of the expansion. In this section, I will briefiy describe effective 
field theories, starting with the pseudoscalars and working up to nucleons 
via the vector particles. I will try to explain some of the physical pictures 
behind the functional dependence from effective field theories, rather than 
provide an exhaustive list of expressions to fit to. 

In Georgi's review |116j of the basic ideas behind effective theory, he 
quotes the essential principles behind effective field theory as 

• A local Lagrangian is used. 

• There are a finite numbers of parameters that describe the interaction 
each of dimension (A; — 4). 

• The coefficients of the interaction terms of dimension {k — 4) is less 
than or of the order of where E < M for some mass independent 
k. 

The above principles allow that physical processes to be calculated to an 
accuracy of i^)'' for a process of energy E. The energy scale M helps to 
organise the power counting. The accuracy of the results can be improved in 
a systematic way if more parameters are included. This is one of the main 
appeals of the formalism. As the energy of the process approaches M, then in 
the strict effective field theory paradigm, a new effective field theory should be 
used with new degrees of freedom ^116j . For example, Fermi's theory of weak 
interactions breaks down as the propagation of the W becomes important. 
For these process the electroweak theory can be used. 

The effective field theory idea can be applied to hadronic physics. In- 
deed historically, the origin of the effective field theory idea was in hadronic 
physics. A Lagrangian is written down in terms of hadron fields, as hadrons 
are the more appropriate degrees of freedom at low energies. The Lagrangian 
is chosen to have the same symmetries as QCD, hence in its domain of validity 
it should give the same physics results as QCD. This is known as Weinberg's 
theorem |12Uj . 

For "large" momentum scales the quark and gluon degrees of freedom will 
become evident, so the theory based on the chiral Lagrangian with hadron 
fields must break down. Unfortunately, it is not clear what the next effective 
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field theory is beyond the hadronic effective Lagrangian, because it is hard 
to compute anything in low energy QCD. 

The lowest order chiral Lagrangian for pseudoscalars is 



V{x) 



(j){x) 
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The expansion parameter |121j for the chiral Lagrangian of pseudoscalar 
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mesons is 



47rF^ ~ IGeV 



(55) 

The power counting theorem of Weinberg [120j guarantees that the tree level 
Lagrangian in equation |S21 will generate the E'^ Lagrangian and non-analytic 
functions. 

The next order terms in the Lagrangian for the pseudoscalars are 
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£4 = ^ OiOi 



(56) 



i=l 



The coefficients are independent of the pion mass. They represent the 
high momentum behaviour. The normalisation used in lattice calculations is 
connected to the one in Gasser and Leutwyler |122j by = 8(47r)^Lj. For 
example, one of the terms in the 0{p^) Lagrangian is 

Oi = [Tt{D^VD^'V^)]^ (57) 

It can be convenient to have quarks with different masses in the sea, 
than those used in the valence correlators. In lattice QCD jargon, this is 
known as "partial quenching" . These partially quenched theories can provide 
information about the real unquenched world |123t I124t I125j . The chiral 
Lagrangian predicts |125t 1126] the mass of the pion to depend on the sea and 
valence quark masses 

]_ ( yujyss - yii) logz/n - y22{yss - ^22) log 1/22 \ 



Mls = yi2{^TxFof[l 



z/22 - yii 

+ ?/i2(2a8 - as) + yssNf{2aQ - a^) 



(5^ 
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GL coeff 


continuum 


nf = 2 [inni 


a5(47rF^) 


0.5 ±0.6 
0.76 ±0.4 


1.22l:l\,{stat)tll 
0.79t%stat)tl\ 



Table 1: where the mass of the sea quark is rus and the masses of the 
valence quarks are mi and m2- Comparison of the results for the Gasser- 
Leutwyler coefficients from lattice QCD (with two flavors) |130j . to non- 
lattice estimates |131j .. 



The function i/ij is related to the quark mass via 



_ i?oK±m,) 

~ (47rFo)2 ^^^^ 



where Fq is the pion decay constant. The leading order term in equation] 
represents the PCAC relation. 

The terms like \ogyu are known as "chiral logs" are caused by one loop 
diagrams. As the chiral logs are independent of the Gasser-Leutwyler coeffi- 
cients (a4, as, ae, as), they are generic predictions of the chiral perturbation 
theory. A major goal of lattice QCD calculations is to detect the presence of 
the chiral logs. This will give confidence that the lattice QCD calculations 
are at masses where chiral perturbation theory is applicable (see j7] for a 
recent review of how close lattice calculations are to this goal). 

There have been some attempts to use equation EHl to determine the 
Gasser-Leutwyler coefficients from lattice QCD [HHl 11211 IHHl IHi • Tabled 
contains some results for the Gasser-Leutwyler coefficients from a two flavour 
lattice QCD calculation at a flxed lattice spacing of around 0.1 fm |130j . 
compared to some non- lattice estimates |131j . 

It is claimed that resonance exchange is largely responsible for the values 
for the values of the Gasser-Leutwyler coefficients jl32[ll3l] . For example, in 
this approach the value of as is related to the properties of the scalar mesons. 



^scalar ^ 8(4vr)2^ (60) 



where Ms is the mass of the non-singlet scalar and the parameters q and c„ 
are related to the decay of the scalar into two mesons. 
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The as and coefficients are important to pin down the Kaplan-Manohar 
ambiguity |133j in the masses chiral Lagrangian at order 0{p^). The latest 
results on estimating Gasser-Leutwyler coefficients using lattice QCD are 
reviewed by Wittig [7|j. 

Most lattice QCD calculations do not calculate at light enough quarks 
so that the log terms in equation EHl are apparent. In practise many groups 
extrapolate the squares of the light pseudoscalar meson as either linear or 
quadratic functions of the light quark masses. For example the CP-PACS 
collaboration |134j used ffi models of the type 

mis = bg^nis + hy^mv + c^^m\ + Cy^niy + Cgymyms (61) 

for the mass of the pseudoscalar {m^'^) where is the mass of the sea quark 
and my is the mass of a valence quark. 

In traditional chiral perturbation theory (see |119j for a review) frame- 
work B is 

Bo = -^{qq) (62) 

It is possible that Bq is quite small |135| 1136] . say of the order of the 
pion decay constant. In this scenario the higher order terms in quark mass 
become important. There is a formalism called generalised chiral pertur- 
bation theory |136j that is general enough to work with a small Bq. The 
chiral condensate has been computed using quenched lattice QCD |137|ll38j . 
The generalised chiral perturbation theory predictions have been compared 
to some lattice data by Ecker |139j . The next generation of experiments may 
be able to provide evidence for the standard picture of a large Bq. 

Morozumi, Sanda, and Soni |14Uj used a linear sigma model to study 
lattice QCD data. Their motivation was that the quark masses on lattice 
QCD calculations may be too large for traditional chiral perturbation to be 
appropriate. 

There is also a version of chiral perturbation theory developed for quenched 
QCD, by Morel |141j . Sharpe fl42j. Bernard and Golterman |143j . Quenched 
QCD is considered as QCD with scalar ghost quarks. The determinant of the 
ghost quarks cancels the determinant of the quarks. The relevant symmetry 
for this theory is f/(3 | 3) x f/(3 | 3). This is a graded symmetry, because it 
mixes fermions and bosons. The chiral Lagrangian is written down in terms 
of a unitary field that transform under the graded symmetry group. 
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In this section, I follow the review by Golterman |144j . The problems 
with quenched QCD can be seen by looking at the Lagrangian for the rj and 
r( (the ghost partner of the rl\ 




\ml{{jif - ijif) + \^\^ - rif + \a\d,{ji - ri)f (63) 



The propagator for the rj' can be derived from equation 1631 



1 



(64) 



The double pole in equation IHH stops the t]' decoupling in quenched chiral 
perturbation theory as /z gets large. This has a dramatic effect on the de- 
pendence of the meson masses on the light quark mass. For example the 
dependence of the square of the pion mass on the light quark mass niq is 



where 6 is defined by 



and A, B, and C are functions of the parameters in the Lagrangian. In 
standard chiral perturbation theory, the S term in equation IHHl is replaced 

2 

by j^p)- As the quark mass is reduced, equation ESI predicts that the pion 
mass will diverge. Values of 5 from 0.05 to 0.30 have been obtained from 
lattice data. It has been suggested |7j that the wide range in 5 might be 
caused by finite volume effects in some calculations. 

The effective Lagrangians, encountered so far, assume that the hadron 
masses are in the continuum limit. In practise most lattice calculation use 
the quark mass dependence at a fixed lattice spacing. The CP-PACS col- 
laboration jnni have compared doing a chiral extrapolation at finite lattice 
spacing and then extrapolating to the continuum, versus taking the contin- 
uum limit, and the doing the chiral extrapolation. In quenched QCD, the 
CP-PACS collaboration [20! found that the two methods only differed at the 
1.5cr level. 



= Amq{l + 5{\og{Bmq)) + Cniq) 



(65) 
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It is very expensive to generate unquenched gauge configurations with 
very different lattice spacings, so it would be very useful to have a formalism 
that allowed chiral extrapolations at a fixed lattice spacing. Rupak and 
Shoresh ^145ij have developed a chiral perturbation theory formalism that 
includes 0(a) lattice spacing errors for the Wilson action. This type of 
formalism had already been used to study the phase of lattice QCD with 
Wilson fermions |146j . An additional parameter p is introduced 

p = 2Woacsw (67) 
where Wq is an unknown parameter and csw is defined in the appendix. This 

2 

makes the expansion a double expansion in -p- and 

= ^tr(5: D^V^D^V) - ^tr((x + p)V + V{x^ + pt)) (68) 

The Lagrangian in equation IHHl is starting to be used to analyse the results 
from lattice QCD calculations [WZ^. 

For vector mesons it less clear how to write down an appropriate La- 
grangian than for the pseudoscalar mesons. There are a variety of different 
Lagrangians for vector mesons |148t I149| 11501 1151j , most of which are equiv- 
alent. A relativistic effective Lagrangian for the vector mesons |15Uj is 

£ = -^Tr(D^K - D,V,f + ^M2Tr(KV^") (69) 

where contains the vector mesons. As discussed earlier in this section, a 
crucial ingredient of the effective field theory formalism is that a power count- 
ing scheme can be set up. The large mass of the vector mesons complicates 
the power counting, hence other formalisms have been developed. Jenkins et 
al. |151j wrote down a heavy meson effective field theory for vector mesons. 

In the heavy meson formalism the velocity v with f ^ = 1 is introduced. 
Only the residual momentum (p) of vector mesons enters the effective theory. 

kv = Myv + p (70) 

In the large A^^^ limit the meson fields live inside A"^. 
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The Lagrangian for the heavy vector mesons (in the large and massless 
limit ) is 

C = -iTTNl{v.V)m - ig2TT{Nl N.jAxv^e'"'^'' (72) 

The connection between the Lagrangian in equation 1^21 and the Lagrangian 
in equation [721 is discussed by the Bijnens et al. |15Uj . 
At one loop the correction to the p mass |151j is 

Am, = -^^L_(gl{'^ml + 2ml + ^^?) + 9i{ml)) (73) 

where gi and g2 are parameters (related to meson decay) in the heavy vector 
Lagrangian and / is the pion decay constant. The next order in the expansion 
of the masses of vector mesons is in the paper by Bijnens et al. |15()j . The 
equivalent expression in quenched and partially quenched chiral perturbation 
theory has been computed by Booth et al. |152j and Chow and Ray |153j . 
The heavy vector formalism suggests that for degenerate unquenched quarks 
the mass of the vector mass should depend on the quark mass nig like 

M^^^ = M^^" + M^^^mg + M^^^mf^ + 0(mj) (74) 

In practise it has been found to be difficult to detect the presence of the 
m^/^ term in equation [73] from recent lattice calculations. The mass of the 
light vector particle from lattice QCD calculations is usually extrapolated to 
the physical point using a function that is linear or quadratic in the quark 
mass. For example the CP-PACS collaboration |134j used the fit model 

j^vec ^ ^vec ^ ^Jec^^ ^ E^^^my + Cj""m| + C^^'^m^ + Cl'^myms (75) 

to extrapolate the mass of the vector meson {M^ec) in terms of the sea 
(mg) and valence {my) quark masses in their unquenched calculations. CP- 
PACS |134j also investigated the inclusion of terms from the one loop calcu- 
lation of the correction to the rho mass in equation [7S1 

There is an added complication for the functional dependence of the mass 
of the p meson on the light quark mass, because in principle the rho can decay 
into two pions (see section [TU)) . This decay threshold complicates the chiral 
extrapolation model. The first person to do an analysis of this problem for 
the p meson in lattice QCD was DeCrand [154 . There was further work done 
by Leinbweber amd Cohen |155j . The effect of decay thresholds on hadron 
masses is also a problem for the quark model |156[ 123,'^ . 
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Figure 6: Self energy contributions to the mass of the rho meson. 

The Adelaide group |157j have studied the issue of the effect of the rho 
decay on the chiral extrapolation model of the rho meson mass in more 
detail. The physical motivation behind Adelaide group's program in the 
extrapolation of hadron masses in the light quark mass has been reviewed 
by Thomas jl58j . 

The mass of the rho is shifted by tt — tt and tc — u intermediate states. 
The effect of the two meson intermediates states can be found by computing 
the Feynman diagrams in figure IHl from an effective field Lagrangian. The 
self energies from the vr — vr {T,^^) and the tt — a; (J^f^^) intermediate states, 
renormalize the mass of the p meson. 

mp = \Jml + S^TT + (76) 



To explain the idea, I will consider the self energy corrections from n — u 
intermediate states in more detail. 

127r2 io w^ikf ^ ' 

where 

w^{kf = e + ml (78) 

and /ip is the physical mass of the p meson. The integral for the S^^ is 
similar, but the algebra is more complicated. 

As the momentum increases the effective field theory description of the 
physics in terms of meson field breaks down. Adelaide prefer to parameterise 
the breakdown of the effective field theory by introducing a form factor Uu)-K{k) 
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at the interaction between the two pseudoscalars mesons and the vector me- 
son. A dipole form factor is used in equation [77| 



--^k) = (|^)^ (79) 

The A^t^ parameter is a (energy) scale associated with the finite extent of 
the hadrons. This is a fit parameter that the Adelaide group |157j determine 
from lattice data. They obtain |157j A^^^ = 630 MeV. 

It is instructive to look at the self energy S^^ with a sharp cut off 

= -'-^(ral arctan( A) + a3 _ ^^2^ (go) 

There is strong dependence on A in equation [HUl The term contains the 
term of chiral perturbation theory. 

In the Adelaide approach |157j the fit model used to extrapolate the mass 
of the rho meson in terms of the quark mass is in equation |^ 

= Co + cml + ^'.i^--,rn.) + KAAn.,rn^) 

2(co + cimi) 

This is a replacement for the chiral extrapolation model in equation 

The Adelaide group |157j note that the coefficient of the term in 
equation ESI is known, hence this is a constraint on the fits from the lattice 
data. However, the actual fits to the mass of the rho particle from lattice 
QCD, do not reproduce the known coefficient in equation [TSl When the 
Adelaide group |157j fit the expression in equation |^ to the mass of the 
rho at the coarse lattice spacings from calculations by UKQCD and 
CP-PACS, the correct coefficient is obtained. 

It is interesting to compare the Adelaide group's |157j approach to the 
chiral extrapolation of the mass of the rho meson to a more "traditional" 
effective field theory calculation. Equation [77| looks similar to a perturbative 
calculation with a cut off. In an effective field theory calculation terms with 
powers of A would be absorbed into the counter terms. In a field theory 
with a cut off, the actual cut off should have no effect on the dynamics 
in the effective theory. A strong dependence on the cut off would signify 
the breakdown of the effective theory. The lecture notes by Lepage discuss 
the connection between a cut off field theory and renormalisation in nuclear 
physics |117j . 
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In his review on effective field theories Georgi quotes Sidney Cole- 
man as asking "what's wrong with form factors?" Georgi's answer is "noth- 
ing". My translation of Georgi's more detailed answer to Coleman's question 
is in the next paragraph. 

Both the Adelaide group's approach and effective field theory agree for 
low momentum scales, hence both formalisms can reproduce the non-analytic 
correction in equation [721 to the mass of the rho, however the two formalisms 
differ in the treatment of the large momentum behaviour. In the effective 
field theory paradigm the large momentum behaviour is parameterised by a 
local Lagrangian with terms that are ordered with a power counting scheme. 
In the Adelaide group's approach the long distance physics is parameterised 
(presumably with Coleman's blessing) by a form factor. In an ideal world, 
the use of an effective field theory is clearly superior to the use of a form 
factors as the accuracy of the approximation is controlled. For the rho and 
nucleon it is not obvious how to set up a power counting scheme in the energy 
(although people are trying). Also an effective field theory based on hadrons 
will no longer describe the physics at large momentum when the quark and 
gluon degrees of freedom become important, hence the use of a form factor 
may be a more pragmatic way to control the long distance behaviour of 
hadronic graphs. It may be easier to introduce decay thresholds in a form 
factor based approach. The hard part of a form factor based approach is 
controlling the errors from the approximate nature of the form factor. The 
Adelaide group jl59j do check the sensitivity of their final results by using 
different form factors. For the light pseudoscalars, the standard effective field 
theory formalism is clearly superior. 

There is a "tradition" of not using "strict" effective field theory tech- 
niques for the p meson. For example DeGrand |154j used a twice subtracted 
dispersion relation to regulate the graph in figure IHl 

Nucleons have also been incorporated into the Chiral Lagrangian ap- 
proach (see |119j for a review). 




(82) 




(83) 
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where rriN and g\ are parameters in the chiral hmit. 



= U (84) 

Ufj, = iv^dfJJv) (85) 

= \[n\d,u] (86) 
D, = d, + T, (87) 

U is the SU (2) matrix containing the pion fields. 

There are a number of comphcations with baryon chiral perturbation 
based on the relativistic action over meson perturbation theory. As discussed 
earlier the power counting is the key principle in effective theory as it allows 
an estimate of the errors from the neglected terms. However, the nucleon 
mass in baryon chiral perturbation theory is the same order as AttFq. This 
complicates the power counting. Also the expansion is linear in the small 
momentum. 

Most modern baryon chiral perturbation theory calculations are done us- 
ing ideas motivated from heavy quark effective field theory |16m I161j . The 
four momentum is factored into a velocity dependent part and a small resid- 
ual momentum part. 

P^l = mv^ + (88) 
The baryon field is split into "large" and "small" fields. 

= e''""-"P,+* (89) 
K, = e^'^^-^P^^^ (90) 

(91) 

where the projection operator is defined by 

P„,,_ = i±|^ (92) 

The leading order Lagrangian for heavy baryon chiral perturbation theory 
(HBChPT) is 

CnN = JIviiv.D + ^-^S.,.u)M (93) 

where 

S^. = '-l^cr^'v. (94) 
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In principle there are corrections to equation IHSl See the review article 
by for a detailed comparison of the relativistic and heavy Lagrangians. 

The final extrapolation formula for the nucleon mass M^v as a function 
of the quark mass rUg is: 



The coefficient C is negative and is a prediction of the formalism. 

The convergence of baryonic chiral perturbation theory is very poor even 
in the continuum. On the lattice the pion masses are even larger, hence 
there are additional concerns about the convergence of the predictions. Using 
chiral perturbation theory, Borasoy and Meissner |162j computed the nucleon 
mass (and other quantities) using heavy baryonic chiral perturbation theory, 
including all quark mass terms up to and including the quadratic order. The 
result for the nucleon mass is in equation IHEl for each order of the quark mass 



where m = 767 MeV. Note this is not a lattice calculation. Although, when 
all the terms are summed up, the correction is small to the nucleon mass, 
there are clearly problems with the convergence of the series. The corrections 
to the masses of the A, E and H baryons were also sizable. 

Donoghue et al. |163j blame the poor convergence of baryoninc chiral 
perturbation theory on the use of dimensional regularisation. Donoghue et 
al. [ 163] argue that the distances below the size of the baryon in the effective 
field theory description breaks down, however the dimensionally regulated 
graphs include all length scales. The incorrect physics from the graphs is 
compensated by the counter terms in the Lagrangian. Unfortunately, this 
"compensation" makes the expansion poorly convergent. Lepage jll7j also 
gives an example (from [164^ ) from nuclear physics where using a cut off gives 
a better representation of the physics than using minimal subtraction with 
dimensional regularisation. 

The graph in equation 1^71 occurs at one loop for the baryon masses. 



= + rrigB + nt'l'^C + .... 



(95) 




(96) 




(97) 



In dimensional regularisation |165j 




(98) 
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The graph in equation |H7| contributes the lowest non-analytic term in the 
octet masses (Mjv oc m^). Donoghue et al. point out that it is a 

bit suspicious that the result in equation EHl is finite when the integral in 
equation inn is cubicly divergent. 

Consider now equation IHZI regulated with a dipole regulator. 

d^k kikj /A(m) 

~*%^T1 — (99) 



(27r)4(A;o-ie)(P-m2 + 2e)'A2-F' ' 247r 



, , 1 .2m + A^ 



In the limit m < A limit 



1 1 

/A(m) — > -A^ - -Am^ + + ... (101) 

Hence for small m the result from dimensional regularisation is reproduced. 

Up to this point the treatment of the integrals looks very similar to the 
approach originally advocated by the Adelaide group |159| I166j . However, 
Donoghue at al. |163j treat A as a cut off. Strong dependence on A is removed 
via renormalisation. 

Mq^ = Mo- dA^ (102) 

where d is a function of the other renormalised parameters in the Lagrangian 
(such as the pion decay constant). In the original work by the Adelaide 
group |159t I166j . the parameter A was a physical number that could be 
extracted from the lattice data. In the formalism of Donoghue at al. |163j . 
the physical results should not depend on A, although a weak dependence 
on A may remain because the calculations are only done to one loop. The 
results for the mass of the nucleon as a function of the order of the expansion 
are 

= 1.143 - 0.237 + 0.034 = 0.940GeV (103) 

with the cut off A = 400 MeV. See |167j for a brief critique of formalism of 
Donoghue at al. |lti3j . 

The Adelaide group |159| 1166] consider the one loop pion self energy to 
the nucleon and delta propagators. The method is essentially the same as 
the one applied to the chiral extrapolation of the rho mass. The Adelaide 
group (159f 166j fit model for the nucleon mass is 

Mn = aN + pNml + aAr7v(m^, A) + aNA{m^, A) (104) 
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where a come from one loop graphs. Equation I1U4I is a three parameter fit 
model: a, /3, and A. 

The Adelaide group have applied the formalism of Donoghue at al. |163j 
to the analysis of lattice QCD data [IHi CHHl from CP-PACS. They have 
compared it with their previous formalism |1681 1168j . 

Lewis and collaborators have studied the lattice regularisation of chiral 
perturbation theory |1691 ll7Uj . This type of calculation can not be used to 
quantify the lattice spacing dependence of lattice results |17Uj . as the lattice 
spacing dependence of the two theories could very different. However, it 
is interesting to explore different regularisation schemes. There is renewed 
interest in the relativistic baryon Lagrangian, because a method |171j . called 
"infrared regularisation" allows a power counting scheme to be introduced 
(see |119j for a review). 

All the above analysis of effective Lagrangians relied on perturbation 
theory to study the theory. Hoch and Horgan |172j used a numerical lattice 
calculation to study the SU{2) x SU{2) non-linear model, for pions and the 
nucleon. The unitary matrix V 

V = j{a.l + niaral^) (105) 
a2 + vr^ = (106) 

S{V) = -V^tr(^r(x)\/(:^ + /i,)t + r(x + /i,)0(x)) (107) 

^ X j = l 

-jFml ^ tT{U{x) + UixY) 

The numerical calculation was done with a small volume 8^, and no finite 
size study was done. The comparison of the lattice results with the pertur- 
bative results is complicated by the effect of the unknown parameters in the 
next order Lagrangian (equation EHlfor example). Hoch and Horgan |172j 
found that the lattice calculation disagreed with the predictions of one loop 
perturbation theory for log divergent quantities. 

A more conservative (some might say cowardly) approach to chiral ex- 
trapolations is to only interpolate the appropriate hadron masses to the mass 
of the strange quark, in an attempt to try to minimise the dependence of any 
results on uncontrolled extrapolation to the light quark masses. One formal- 
ism 11735 for doing this is called the "method of planes". Similar methods 
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have been used by other groups (see for example |174| I175j ). Obviously, this 
type of technique is not useful to get the nucleon mass. In unquenched QCD, 
the sea quark masses should be extrapolated to their physical values, so there 
is no way to avoid a chiral extrapolation even for heavy hadrons. 

It is traditional to plot the hadron masses before any chiral extrapolations 
have been done, so as not to contaminate the raw lattice data from the 
computer with any theory. In an "Edinburgh plot" the ratio of the nucleon 
to rho mass is plotted against the pion to rho mass |176j . If there were no 
systematic errors, such as lattice spacing dependence, then the data should 
fall on a universal curve. It is also common to use 1/ro [12] (see section 
13. 2j) as a replacement for the mass of the p. There is also an APE plot that 
plots the ratio of the nucleon to rho mass against the square of the pion to 
rho mass |177j . This parameterisation is meant to have a smoother mass 
dependence than the Edinburgh plot. 

3.4 FINITE SIZE EFFECTS. 

The physical size of the lattice represents an obvious and an important sys- 
tematic error. One simple way to estimate the size of a hadron is to consider 
its charge radius. For example, the proton's charge radius is quoted as 0.870 
(8) fm in the particle data table jH]. The sizes of all recent unquenched lattice 
QCD calculations are all above 1.6 fm (see sectional). The fact that all the 
physical lattice sizes were much bigger than the charge radius does not rule 
out finite size effects. In this section I will discuss some of the mechanisms 
thought to be behind finite size effects in lattice data. 

A nice physical explanation of the origin of finite size effects for hadron 
masses has been presented by Fukugita et al. |178j . Consider a hadron in 
a box with length L and periodic boundary conditions. Most lattice QCD 
calculations use periodic boundary conditions in space. The path integral for- 
malism requires that the boundary conditions in time are anti-periodic |179j . 
although many groups use periodic boundary condition in time as well as 
space. The finite size of the box will mean that a hadron will interact with 
periodic images a distance L away. The origin of finite size effects is closely 
related to nuclear forces. The MILC collaboration have described a qualita- 
tive model for finite size effects based on nuclear density [121 • The nucleon 
is considered as part of nucleus comprising of the periodic images of the 
nucleon. 
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The self energy of the hadron (SE) will be 



6E = Y,V{nL) (108) 

n 

where V{x) is the potential between two hadrons a distance x apart. Various 
approximations to V{x) give different functional forms for the dependence of 
the hadron mass on the volume. 

In the large L limit the potential is approximated by one particle exchange 
(e~™'^/r). The interaction energy goes like V{0)+6V{L), so the dependence of 
the hadron mass will be ~ g-"ijri This argument can be made rigorous jlSUI 

CHU. 

It is useful to consider equation llOSl in momentum space using the Poisson 
resummation formulae. 

5E=^j:nn.Y^ (109) 

A more general expression for the potential between two hadrons can be 
derived if the spatial size of the hadron is modelled with a form factor F(k). 

= (110) 

Momentum is quantised on the lattice in quanta of In physical units the 
value of a quantum of momentum is around 1 GeV, hence the k = term in 
equation 11091 should dominate the sum. Therefore this model predicts that 
the masses should depend on the box size like: 

M = M^ + cL-^ (111) 

This model is physically plausible but not a rigorous consequence of QCD. 
Fukugita et al. |178j noted their data could also be fit to the functional form. 

m^ = ml + c^ (112) 

rather then equation lllll Unfortunately only the theory for the regime of 
point particle interacting at large distances is really rigorous |180t I18lj .. 

Chiral perturbation theory can be used to compute finite volume cor- 
rections. For example the ALPHA/UKQCD collaboration have used |182j 
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a chiral perturbation theory calculation by [183 Gasser and Leutwyler to 
estimate the dependence of the pion mass mps{L) on the box size L. 



mps{L) 1 m 



l + ^'^g{mpsL) (113) 
mps[oo) Nf Fps 

= ^ r ^ e-^^^-V(-) (114) 

Garden et al. |182j used equations II 131 and II 141 to show that the error in 
mps was 0.1 % for mpsL > 4.3. This an example of the rule thumb that finite 
size effects in hadron masses become a concern for irtpsL < 4. Golangelo et 
al. |184j are extending equation 11131 to the next order. Ali Khan et al. jl85j 
have started to use chiral perturbation theory to study the finite size effects 
on the nucleon. 

It would be useful if some insight could be gained from finite size effects in 
quenched QGD that could be applied to unquenched QGD where the volumes 
are smaller. Unfortunately, there are theoretical arguments |186j that show 
that finite size effects in the unquenched QGD will be larger than in the 
quenched QGD. A propagator for a meson can be formally be written in 
terms of gauge invariant paths. Gonceptually this can be understood using 
the hopping parameter expansion. The quark propagators are expanded in 
terms of the k parameter. The hopping parameter expansion |187j was used 
in early numerical lattice QGD calculations, but was found not to be very 
convergent for light quarks. 

c{t) = E ^'S\W{C)) + E ^'P^UP{C')) (115) 
c c 

where W{C) are closed Wilson loops inside the lattice of length L{C). P{C') 
are Polykov lines that wrap around the lattice in the space direction. The 
aval parameter is 1 for periodic boundary conditions and —1 for anti-periodic 
boundary conditions for Polykov lines that wrap around the lattice an odd 
number of times. 

In quenched QGD {P{C')) = because of a symmetry of the pure 
gauge action, so only the first term in equation 11151 contributes to the cor- 
relator c{t). The centre of SU(3), the elements that commute with all the 
elements of SU(3), is the Z{3) group [T7j. The Wilson gauge action is invari- 
ant under the gauge links being multiplied by a member of the centre of the 
group on each time plane. In unquenched QGD, the Z-^ symmetry is broken 
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by the quark action, so both Wilson loop and Polyakov lines contribute to 
the correlators. There is no clear connection between the arguments based 
on equation 11151 and the nuclear force and chiral perturbation formalism for 
finite size effects. 

In their large scale quenched QCD spectroscopy calculations the CP- 
PACS collaboration kept the box to be 3 fm so did not apply any 
corrections for the finite size of the lattice. Using the finite size estimate from 
previous calculations CP-PACS [^2] estimated that the finite size effects were 
of the order 0.6 %. The GFll group ^ did a finite size study at a coarse 
lattice spacing (a = 0.17 fm). The results for ratios of hadron masses at finer 
lattice spacings were then extrapolated to the infinite volume limit using 
the estimated mass shift at the coarse lattice spacing. Gottlieb j^Zj gives a 
detailed critique of the method used by the GFll group to extrapolate their 
mass ratios to the infinite volume. 

If a formalism could be developed that would predict the dependence of 
hadron masses on the box size, then this would help make the calculations 
cheaper (see equation HUj). The additional savings in computer time could 
be spent on reducing the size of the light quark masses used in the calcula- 
tions |188j . Gottlieb suggested, the only way to control finite size effects 
is to keep increasing the box size until the masses no longer change. 

The finite box size is not always a bad thing. The size of the box has 
been used to advantage in lattice QCD calculations (see |189j for a review). 
A definition of the coupling is chosen that is proportional to 1/L (where 
L is the length of one side of the lattice). A recursive scheme is setup that 
studies the change in coupling as the length of the lattice side L is halved. The 
Femtouniverse, introduced by Bjorken |190j . is a useful regime to study QCD 
in. There is a chiral perturbative expansion based on the limit 171^^ i; 
(see for a modern application). See Van Baal |191j for a discussion of 
the usefulness of the finite volume on QCD. 

It seems possible to run calculations in a big enough box to do realistic 
calculations. So the prospects are good for a first principles lattice calculation 
of hadron masses without resort to approximations in simulations of that 
kind required in condensed matter systems in simulations of macroscopic 
size systems. 
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4 AN ANALYSIS Of SOME LATTICE DATA 



To consolidate the previous material, I will work over a simple analysis of 
some lattice QCD data. I think it is helpful to understand the steps in lattice 
QCD calculation, if the ideal case where QCD could be solved analytically is 
considered first. All the masses of the hadrons would be known as a function 
of the parameters of QCD: quark masses (mj) and coupling g. 

Mh = fuimu, ma, m^, nic, rrib, g) (116) 

As the masses and coupling are not determined by QCD, the equation for 
all the hadrons H would have to be solved to get the parameters. The 
solution would be checked for consistency that a single set of parameters 
could reproduce the entire hadron spectrum. For a calculation of a form 
factor, the master function fu would also depend on the momentum. 

The formalism of lattice QCD is in some sense is a discrete approximation 
to the function fn- The result from calculation would be a table of numbers: 

Mh = fu^mu, ma, rus, rUc, rUb, g, L, a) (117) 

An individual lattice calculation would also depend on the lattice spacing 
and lattice volume. Lattice calculations have to be done at a number of 
different lattice spacings and volumes to extrapolate the dependence of /|^** 
on L and a. 

By doing calculations with a number of different parameters the results 
can be combined to produce physical results in much the same way that could 
be done if the exact solution was known. In particular the lattice spacing 
and lattice volume have be extrapolated away to get access to the function 
fn. 

To understand the procedure in more detail, I will work through a naive 
analysis of some lattice QCD data from the UKQCD collaboration |192j . 
Table |21 contains the results for the mass of the vr and p particles in lattice 
units from a quenched QCD jl92j . The lattice volume was 24^ 48, [3 = 
6.2, and the ensemble size was 216. The clover action using the ALPHA 
coefficients was used. To use the data in table El the secret language of the 
lattice QCD cabal must be converted to the working jargon of the continuum 
physicist. 

The ensemble size of 216 means that 216 snapshots of the vacuum (the 
value of in equation ITS|) were used to compute estimates of the p and 
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arrip 


0.13460 
0.13510 
0.13530 


0.2149il^ 
0.1836t?| 


0.3887li 
0.353lii 
0.34141^^ 



Table 2: Lattice QCD data for the masses of the pion and rho from lattice 
QCD calculations from UKQCD [TH^ 



TT correlators. A supercomputer was used to compute the correlators for 
each gauge configuration from quark propagators (see equation HH)). The 
masses were calculated by fitting the correlator to a fit model of the form 
in equation 1221 using a minimiser such as MINUIT |193j . The error 
bars in the table in equation |21 come from a statistical procedure called the 
"bootstrap" method [M]. 

A table of numbers of hadron masses (or even a graph) is not too useful. 
A better way to encapsulate the hadron masses as a function of quark mass 
is to use them to tune effective Lagrangians, as discussed in section 13.31 To 
plot the data in table |21 in a more physical form, I convert from the n value 
to the quark mass 

m, = i(i-J-) (118) 

There are additional 0{a) corrections to equation lll8l |185] . The Hcru param- 
eter is required because clover fermions break chiral symmetry. The value 
of Kcrit is the chosen to give a zero pion mass. Equation 11181 is the basis of 
the computation of the masses of quarks from lattice QCD. However, pertur- 
bative factors are required to convert the quark mass to a standard scheme 
and scale. This perturbative "matching" can be involved, so the value of the 
quark mass is rarely used to indicate how light a lattice calculation is. 
The simplest thing to do is to use a fit model in equation 11191 

i^ps = Spsruq (119) 

There are classes of fermion actions (see appendix such as staggered or 
Ginsparg- Wilson actions, that do not have an additive renormalisation. 
A simple fit to the data in table |21 with the fit model in equations 11191 
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Figure 7: Rho and pion masses in lattice units as a function of inverse kappa 
value. The data set had (3 = 6.2, andavolumeof 24^48. The fermion operator 
was the non-perturbatively improved clover operator. 



and lllSl gives Kcrit = 0.135828(6), to be compared to Kcru = 0.135818l|7 from 
the explicit analysis from UKQCD that included a 0{a) correction term. 

To simplify the analysis, I will assume that the physical mass of the light 
quark is zero. I fit the data for the mass of the rho in table El to the model 
in equation I12()[ 

my = Ay + Syrriq (120) 

The result for the Ay parameter is 0.304(2). If the mass of the light quark 
is assumed to be zero, then anip = 0.304, hence the lattice spacing is 2530 
MeV (using nip = 770 MeV). This can be compared against the = 2963 
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MeV from tq [Ml- 

As the masses of quarks get lighter more sophisticated fit models based 
on the ideas in section 13.31 can be used. Although the basic ideas behind 
the analysis outlined in this section are correct, there are many improve- 
ments that can be made, particularly if the rho and pion correlators for each 
configuration are available. 

5 PARAMETER VALUES OF LATTICE QCD 
CALCULATIONS 

The results from unquenched lattice QCD calculations, with the lattice spac- 
ing and finite size effects accounted for, are the results from QCD at the 
physical parameters of the calculation. Hence a key issue in unquenched 
calculations is how close the parameters are to the physical parameters. For 
example, as discussed in section 1^51 ideally the masses of the quarks must be 
light enough to match the lattice results to chiral perturbation theory. The 
parameters used in a lattice calculation are usually dictated by the amount of 
computer time available, or what gauge configurations are publicly available. 

In this section, I will describe the current state of the art in the parameters 
used in lattice QCD calculations of the hadron spectrum. It is not entirely 
obvious which parameters to use to characterise a lattice calculation. The 
obvious choice of using quark masses is complicated by the need for renor- 
malisation and running. To show how light the quarks are in a calculation, I 
plot the ratio of the pseudoscalar mass to the vector mass for as a function of 
lattice spacing. The danger of this type of plot is that it says nothing about 
finite size effects. I usually just show the ratio for the lightest quark mass, 
as this is the most computationally most expensive point. The error bars on 
the ratio gives some indication on the statistical sample size. I have always 
used a lattice spacing defined by tq = 0.49 fm (see section ing . 

In figure IHl I plot the smallest pseudoscalar mass to vector mass ratio as a 
function of lattice spacing for some recent large quenched calculations from 
the following collaborations: GFll ij, CP-PACS [1201, UKQCD [TT?^ 
BGR [197J, and ALPHA jl98j . The improvements in parameters between 
the GFll calculation [SlIU, at the start of 1990's and the CP-PACS calcula- 
tion |196j at the end of 1990's can clearly be seen in the figure. 

Both the main benchmark calculations by CP-PACS |19tij and GFll jSlH] 
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Figure 8: Smallest ratio of pseudoscalar mass to vector mass as a function 
of lattice spacing for a number of quenched QCD calculations. The labels 
have the collaboration name followed by the type of fermion operator. 

used both the Wilson fermion and gauge actions. As I discuss in section El 
lattice QCD actions with better properties have been developed. For ex- 
ample, the point in figure |H1 from UKQCD used the clover fermion action. 
The clover action is designed to have smaller lattice spacing errors compared 
to the Wilson action. Unfortunately, the calculations by UKQCD |192j and 
QCDSF 199J that used the non-perturbative improved clover action reported 
problems with exceptional configurations, hence for those lattice spacings the 
light quark masses can not be reduced further at these lattice spacings. There 
is a formalism called twisted mass QCD |198j that is a natural extension of the 
clover action, that can be used to explore the light quark mass regime |2()()j . 
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Collaboration 




a fm 


L fm 


Mps 
Mv 


mps MeV 


MILC ^ 


2+1 


0.09 


2.5 


0.4 


340 


CP-PACS nn 


2 


0.11 


2.5 


0.6 


900 


UKQCD [73 


2 


0.1 


1.6 


0.58 


600 


SESAM pnHll^n4 


2 


0.074 


1.8 


0.57 


530 



Table 3: Typical parameters in recent unquenched lattice QCD calculations. 

To reach lighter quark masses, it seems likely, that new fermion actions 
such as the overlap- Dirac, Domain wall, fixed point [201 , and twisted mass 
QCD |198j will be required to study the light quark mass region of quenched 
QCD, and hence improve on the results of the CP-PACS calculation |196j . 

The first results from these types of calculations are shown in the figure |H1 
Unfortunately, the fermion operators that obey the Ginsparg- Wilson relation 
are computationally more expensive than Wilson like actions. At the moment 
the error bars are too large to be competitive with those from actions that 
use Wilson type fermions. However the control of the systematic errors in the 
calculations that use Ginsparg- Wilson operators are rapidly improving |197j . 

TableElshows the parameters of some recent large scale unquenched calcu- 
lations. I have also included the mass of the lightest pion in the calculations. 
Although the choice of the lattice spacing in a specific calculation can have 
large uncertainties, I feel that the mass of the pion in physical units gives 
a more immediate measure of the "lightness" of the quarks in a calculation. 
For the table I just used the lattice spacing quoted by the collaboration. 
The review article by Kaneko 1^1^ contains a more thorough survey of the 
parameters of some recent unquenched calculations. 

In figure El the ratio of the of masses of the lightest pseudoscalar to vector 
for some recent full QCD calculations is plotted. The effect of the cost for- 
mula for full QCD^nican be clearly seen, as the finest lattice spacing used is 
around 0.1 fm. Only the CP-PACS calculation had a number of different lat- 
tice spacings. The non-perturbative improved clover action was constrained 
to lie around 0.1 fm, because the computation of the clover coefficient was 
too expensive at coarser lattice spacings 205^ . 

The most interesting point is from the calculations done by the MILC 
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Operator 


jPC 


Lightest particle 




0-+ 


TT 


V'l7475V'2 


0-+ 


TT 


^^17*^2 


1— 


P 


^l747iV'2 


1— 


P 


^i7i7jV^2 


1+- 


hi 


^17*75^2 


1++ 


ai 


^1^2 


0++ 


ao 



Table 4: Interpolating operators for light mesons. The 1 and 2 subscripts 
label flavour and show that the mesons are non-singlet. 

collaboration (HH]. The MILC collaboration are currently running with a 
lattice spacing of 0.09 fm. 

6 THE MASSES OF LIGHT MESONS 

Light mesons have a number of important uses in lattice QCD calculations. 
In calculations that use Wilson like fermions, the mass of pion is used to 
calculate the additive mass renormalisation {kctu in equation II 18|1 . The p is 
sometimes used to set the lattice spacing. The mass of one of the mesons: 
kaon, i^*, or is used to calculate the strange quark mass. After the light 
quark masses are calculated any remaining masses are used as a consistency 
check. The masses of the quarks are used in any further calculations, such 
as the computation of matrix elements. 

The interpolating operators for mesons to be used in equation ^2 are in 
table 0] The J^^ quantum numbers of the meson operators can be derived 
using the standard representation of the parity P and charge conjugation 
operators from the Dirac theory. The meson interpolating operators are 
usually extended in space using one of the prescriptions in section |21 Most 
calculations concentrate on the S-wave mesons as the signal to noise ratio is 
better for these mesons, than for P-wave states. I discuss P-wave states in 
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Figure 9: Smallest ratio of pseudoscalar mass to vector mass as a function 
of lattice spacing for a number of full QCD calculations. The data is from: 
CP-PACS [202J (diamonds), UKQCD [73j (bursts), MILC [90j (octagon), and 
SESAM ^ (squares). 
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section 16.11 

The QCD field strength tensor has also been used with the fermion bilin- 
ears in tabled] The QCD field strength tensor (Fjk) has specific J^*" quantum 
numbers that can be used to obtain fermion bilinear operators with different 
J^*" quantum numbers. For example the MILC collaboration |2U6j used an 
interpolating operator of the form: 

€ijk^i7iFjkip2 (121) 

with J^^ = 0-+. The MILC collaboration obtained the same mass for 
the pion using the operator in equation 11211 as the pion operator in table 0] 
The operator in equation 11211 is more "gluey" than a ip'-f^ip operator, so 
might be expected to couple to hybrid (quark-antiquark pair with excited 
glue) pions. However, the pion is still the hghtest state that couples to the 
operator in equation 11211 hence is the state extracted from the fits. 

In figure ^1 I plot the mass of the vr and p mesons as a function of the 
physical box size for Wilson fermions at a fixed lattice spacing {(3 = 6.0, a ~ 
0.1/m). The data seems to show that a linear box size of 3 fm is big enough 
for "small" finite size errors. The errors on the masses of the mesons with 
the lighter quark masses are too big too draw any conclusions. 

To get a more quantitative estimate of the volume dependence UKQCD |192j 
studied finite size effects in quenched QCD at f3 = 6.0 with non-perturbatively 
improved clover fermions. Two volumes were used with sides 1.5 fm and 3 
fm. For the pseudoscalar channel there was a 2a difference between the mass 
on the two volumes. There were no statistically significant difference between 
the masses for the vector particle between the two lattice volumes. 

To give some idea of the size of lattice spacing errors in the masses of 
the light mesons, I plot the mass of the K* meson as a function of the 
lattice spacing in figure ITT] from the CP-PACS collaboration 96j. The figure 
shows the difference between using the kaon or (p meson to determine the 
strange quark mass. The clear difference between the mass of the K* when 
the strange quark mass is determined in two different ways is caused by the 
quenched QCD not being the theory of nature. The results for meson masses 
in the continuum limit are in table |96j . 

The physical summary of the results in table El is that the hyperfine 
splitting is too low from quenched lattice QCD. To isolate the reduction 
of the hyperfine splitting, Michael and Lacock jl74| 1175] introduced the J 



50 



2 







O 


k=0.1530, pion 






□ 


k=0.1550 pion 






o 


k=0.1563, pion 


1.5 




o 


k=0.1530,rho 






A 


k=0.1550,rho 






V 


k=0.1563,iho 



0.5 



B B 



ii 



1.5 



2 

Lfm 



2.5 



W 
B 



Figure 10: Dependence of the mass of the vr and p mesons on the box size 
in quenched QCD for Wilson fermions at /? = 6.0. The data was taken from 
the compendium of World data in 

parameter. 

J = ^j^2 \mv=i.sMp (122) 

where My and Mp are the vector and pseuodscalar masses respectively. The 
condition My = l.SMp corresponds to the experimental ratio of K* and K 
masses. This mass ratio was chosen so that an extrapolation to quark masses 
below strange was not required. Some theoretical problems with the value 
of J defined at a light quark reference point are discussed by Leinweber et 
al. [TH7] . 
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Figure 11: Lattice spacing dependence of the K* meson in quenched QCD 
from the CP-PACS collaboration 



The J parameter has been chosen to be independent of the lattice spacing 
and a explicit definition of the quark mass. One experimental estimate for J 
is obtained from 

J = mK*^ ^ = 0.48 123 

Including the uncertainty from using (p — K* difference rather than the K* — p 



mass difference, Michael and Lacock ll_74| estimate the experimental value 
of J to be 0.48(2). 

In quenched QCD, Michael and Lacock |174j obtained J = 0.37(2) (4) in 
disagreement with the experimental estimate. Using their bigger data set, 
the CP-PACS collaboration jHEl obtained J = 0.346(23) from a continuum 
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Mass 


Result {rriK ) GeV 


Result (m^ ) GeV 


Expt GeV 


rriK 




0.546(06) 


0.496 


rriK* 


0.846(07) 


0.891(05) 


0.892 




0.970(06) 




1.020 



Table 5: Masses of light S-wave mesons in quenched QCD from CP- 
PACS jnni- The different analyses depend on which meson is used to de- 
termine the strange quark mass. 

extrapolation in quenched QCD. 

The clear disagreement of the J parameter from quenched QCD with 
experiment, in principle makes J a good quantity to measure the effect of 
the sea quarks. The actual definition of J in equation 11221 for unquenched 
QCD is not trivial because, as discussed in section 13.21 the lattice spacing 
does depend on the mass of the sea quarks. The issues in defining J for 
unquenched QCD are discussed by the UKQCD collaboration [73] . 

Figure IT^ compares the value of J from the unquenched lattice QCD from 
UKQCD [IS] and MILC Both the calculations done by the UKQCD and 
MILC collaborations were at a fixed lattice spacing. The calculation by MILC 
was done at a lighter sea quark mass then UKQCD, that presumably explains 
why the value of J from MILC agrees better with the experimental value. 
The MILC result needs to be confirmed by a study of the lattice spacing 
dependence. Kaneko has recently reviewed the status of calculations of 
J from unquenched calculations. 

The effect of unquenching on the hyperfine splitting in light mesons has 
an important effect on the value of the strange quark mass extracted from 
lattice data. For example, the CP-PACS collaboration |202j have used nf = 2 
lattice QCD to calculate the strange quark mass. CP-PACS's results are in 
table El The results show a sizable reduction in the mass of the strange quark 
between quenched and two flavour QCD. Lubicz [20 7j and Wittig |7i review 
the determination of the masses of quarks from lattice QCD. 
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Figure 12: Comparison of the J parameter from UKQCD |73j (diamonds) 
and MILC (HHl (bursts). The experimental point is the square. 
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Table 6: Mass of the strange quark from CP-PACS |2()2j in the MS scheme 
at a scale of 2 GeV. 

6.1 P-WAVE MESONS AND HIGHER EXCITATIONS 

The calculation of the masses of ground state mesons and baryons is essen- 
tially done to validate lattice QCD methods and to compute quark masses 
and the strong coupling. For excited S-wave and P-wave states, the issues 
are to determine the quark and glue content of the state (as well as to check 
that the masses agree with experiment). The signal to noise ratio is worse for 
P-wave mesons than for S-wave mesons [72], so the calculations are harder. 
These type of calculations tend to be more exploratory, so usually no attempt 
is made to take the continuum or infinite volume limit. As the underlying 
lattice techniques mature, there will be attempts quantify all the systematic 
errors. 

The interpolating operators in table E] can be used to create the P-wave 
ao, &i, and ai states. However, in the quark model these states have zero wave 
function at the origin. So non-local interpolating operators with a node at 
the origin have been tried '2081 120911210) . Reasonable results can be obtained 
using the operators in table 0] Spin 2 states are not accessible from local 
interpolating operators, so must be created using non-local operators. 

A more general approach |211j is to consider a generic non-local meson 
operator at a specific time t (as is required for an interpolating operator for 
a time correlation function). 

0(r)=^(x,t)r n Ui){x + r,t) (124) 

x,x+r_ 

where the quark and anti-quark are separated by a distance r and F ia an ar- 
bitrary gamma matrix. The IIx.x+t- ^ is a set of gauge links that connects the 
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J=3 


A2,TuT2 


J=4 


Ai,E , Ti,T2 



Table 7: Representation of SU{2) in terms of representations of the cubic 
group. 

quark to the anti-quark in a gauge invariant way. In general Hx.x+r U is not 
unique and will effect the transformation properties of 0{r) under the cubic 
group. Operators are designed to transform under specific representations of 
the cubic group. 

The connection between the representations of the cubic group and the 
SU (2) rotation group are in table [7| The dimensions of the representations: 
Ai, A2, E, Ti, T2 are 1,1,2,3, and 3 respectively. Hence the dimensions of the 
representations match between the cubic group and SU (2) rotation group 
with representation of 2 J + 1 in table 

If the gauge configurations are fixed to Coulomb gauge then non-local 
interpolating operators based on spherical harmonics can be used. This ap- 
proach was studied by DeGrand and Hecht |212[ 17^ . As hadron masses are 
gauge invariant quantities, gauge fixing at an intermediate stage should not 
effect the final results. Meyer and Teper discuss how to construct higher spin 
glueball operators |213j . 

There are number of interesting puzzles with the phenomenology of "P- 
wave" mesons. For example, there are speculations that the ao(980) particle 
is potentially not well described by a qq state in the quark model. Its mass 
is very close to the threshold for two kaon decay. There are many models 
that treat this state as a KK molecules [214j or Tjqqq state. As lattice QCD 
calculations use the non-singlet qq operator to create this state, it may not 
couple strongly to a molecular qqqq state. Therefore I would expect that the 
lightest state in the qq channel to be the a0(1450). However, this speculation 
should, and will be |215j tested in unquenched lattice calculations. 
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As discussed in section I3.H the interpretation of the Qq state is com- 
phcated by a quenched chiral artifact [Sn] in quenched QCD. The largest 
systematic study of the oq in quenched QCD was done by Lee and Wein- 
garten j86j. Unfortunately, they could only find the mass of the oq mass at 
the mass of the strange quark, possibly because of problems with the artifact 
in this channel. Alford and Jaffe |lU9j review some of the earlier results for 
the mass of the Oq from quenched QCD. In a lattice calculation that used 
domain wall fermions and modelled the quenched chiral artifact, Prelovsek 
and Orginos |216j obtained the lightest state in the oq channel to be 1.04(7) 
GeV. At one coarse lattice spacing, Bardeen et al. jHSj obtained a value of 
1.34(9) GeV for the lightest state in the channel, also with an analysis 
that was aware of the quenched chiral artifact. The existing quenched lattice 
QCD data does not determine the mass of the lightest non-singlet O"*""*" state 

In a two flavour unquenched calculation the UKQCD collaboration |215j 
quote the prehminary result for the mass of the ao to be 1.0(2) GeV at one 
lattice spacing. The MILC collaboration jSH] have computed the Oq mass with 
2+1 flavours of sea quarks. The lightest mass of the Oq in MILC's calculation 
is 0.81 GeV j^. The MILC collaboration |HH] also claimed to see evidence 
for the open decay channels of the oq state. Hence, more work is required 
in unquenched QCD to determine the mass of the lightest non- singlet O"*""*" 
state. 

The only values of the masses of the Oi and hi mesons from unquenched 
QCD, I could find were from the MILC collaboration jS^. MILC obtained 
masses of the ai and hi of 1.23(2) and 1.30(2) GeV compared to the exper- 
imental values of 1.23(4) GeV and 1.2295(32) GeV respectively 0. There 
doesn't seem to be any interesting experimental issues with the ai and hi 
mesons. As noted by MILC jHSjj according to the quark model the masses 
of the ai and hi mesons should be similar to the mass of the Oq. Hence, 
any large splitting between the masses of oq and ai and hi is indicative of 
dynamics beyond the simple quark model. 

The most surprising aspect of hadron spectroscopy are the existence of 
Regge trajectories. Empirically, the square of the mass of a hadron with 
mass m{l) is linearly related to the spin by: 

l = a'M'^{l) + a{Q) (125) 

Although the Regge trajectories can be explained by models (see |217j for a 
review), in particular the string model, the existence of Regge trajectories 
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has not been shown to be a rigorous consequence of QCD. Equation 11251 is 
a very useful tool in classifying baryon states. Historically some of the spin 
assignments of baryons were guessed from equation 11251 [218"' . Equation 11251 
is also useful in scattering experiments (see [219ij for a discussion). 

The simple relation in equation 11251 has recently been challenged by a 
number of authors. The improved precision of the experimental data on 
hadron masses has allowed fits to the spectrum that seem to show some 
nonlinearity in the relation between the square of the meson masses and their 
spin \220\ . Brisudova et al. |221j discuss various hadronic models that can 
reproduce linear and nonlinear Regge trajectories. The nonlinearities in the 
Regge trajectories were related to flux tube breaking |217j . hence the study 
of Regge trajectories is complimentary to the study of string breaking |222j . 
Brisudova et al. |217j make the prediction that are no light quark quarkonia 
beyond 3.2 GeV due to the termination of the Regge trajectories. If this 
conjecture is correct, then this might simply searches for hybrid and glueballs 
in the region from 3 to 5 GeV. 

The problem for lattice QCD calculations is that it is hard to construct 
interpolating operators with spin higher than 2 on the lattice. As higher 
orbital states have larger masses than ground state mesons, they have a worse 
signal to noise ratio, so it is hard to extract a signal from the calculations. 

There have been some calculations of spin 3 masses from lattice QCD 
calculations. The UKQCD collaboration |211j studied spin 2 and spin 3 
states in a quenched lattice QCD calculation. In figure IT^ I plot their results 
in a Chew-Frautschi plot |223j with experimental data from the particle data 
table. 

There have been a few attempts to study the spectroscopy of excited 
mesons using lattice QCD. As discussed in section |21 the computation of 
the masses of excited states requires a multi-exponential fit to the lattice 
correlators that can be unstable. 

The CP-PACS collaboration [7^] used the maximum entropy method 
(briefly revived in section 12)) to study the excited spectrum of the rho and 
pion mesons. The calculations were done in quenched QCD, using the same 
data set that was used for their calculation of the ground state masses |196j . 

A spectral density from CP-PACS [7^ is in figure El The masses of 
the states are obtained from the peaks in figure Figure El from CP- 
PACS jTH] shows the masses of the excited p meson as a function of lattice 
spacing. The diverging graph in figure [T31 is thought to be a bound state of 
fermion doublers. 
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Figure 13: Square of the mass of the meson mass versus its spin /. The 
octagons are the experimental numbers and the bursts are the lattice QCD 
data from UKQCD |211j . The quarks have the mass of the strange quark. 
There are only error bars on the lattice data. 

The final results from CP-PACS jT^] were that the first excited state of 
the pion had a mass of 660(590) MeV and the mass of the first excited rho 
meon was 1540 (570) MeV from quenched QCD. The errors also include an 
estimate of the error from taking the continuum limit. 

Experimentally |6j, the first excited pion is the 7r(1300) with the mass 
of 1300 ± 100 MeV. The excited states of the p meson are more interesting. 
There are two states: p(1459) with a mass of 1465 ±25 and the p(1700) with 
a mass of 1700 ± 20 MeV. Donnachie |224j reviews the evidence for a hybrid 
state in the 1 channel. 
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Figure 14: Spectral function (from MEM) from the CP-PACS collaboration 
for the p correlators. The lattice spacing is = 2.58 GeV. The Kll, K42, 
and K44 keys are the plots for mesons with different quark masses. 

The errors on the excited states from CP-PACS are really too big for a 
meaningful comparison with experiment. The results from CP-PACS [71|] 
need confirmation from other, less Bayesian based fitting methods, such as 
variational methods. 



7 THE MASSES OF LIGHT BARYONS 

An important, but perhaps slightly boring, goal of lattice QCD is to compute 
the mass of the nucleon with reliable errors from first principles. The nucleon 
is the most important hadron for the real world of the general public, but the 
nucleon's role in the esoteric domain of particle physics is as "background" 
to the more interesting stuff. In this section, I will discuss the highlights of 
some recent large scale quenched and unquenched QCD calculations. As the 
aim of lattice QCD calculations for these baryons is to validate lattice QCD 
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Figure 15: Continuum extrapolation of the masses in the vector channel 
from CP-PACS. 

before the spectroscopy of more interesting hadrons is attempted, the focus 
will be on the error analysis. 

In principle the mass of the nucleon should be an ideal quantity to com- 
pute on the lattice, because it is stable within QCD, hence there are no 
concerns with the formalism due to decay widths. The main complication 
with getting an accurate value of the nucleon mass is the large chiral extrap- 
olation required (see section I3.3p caused by the large quark masses used in 
the calculations. In the majority of lattice QCD calculations, electromag- 
netism is ignored, so I refer to the generic nucleon rather than the neutron 
or proton. 

The "standard" interpolating operators for the nucleon 1225j 

Ni"^^ = eijk{uJC-f5dj)uk 

where u and d are operators that create the up and down quark respectively 
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and C is the charge conjugation matrix. In the non-relativistic hmit (keeping 
the upper components) the N2 operator vanishes, so it doesn't couple strongly 
to the nucleon. Empirically the N2 operator has been found useful for N* 
states. Leinweber [226 discusses the connection between the interpolating 
operators for the nucleon used on the lattice and those used in QCD sun 
rules. 

The nucleon correlator is constructed, in the standard way, by creating a 
nucleon at the origin, and then destroying it at time t later. 

C±{t) = ^(0 I iV(x,t)(l ±74)iV(0,0) I 0) 

x_ 

The specific representation of equation l23|l is slightly more subtle. 

where m_|_ and m_ are the masses of the lightest positive and negative parity 
states. 

To illustrate the finite size effects on the mass of the nucleon in quenched 
QCD, I plot the nucleon mass as a function of box size for Wilson fermions 
at /? = 6.0 in figure UHl This data set was chosen because there was data 
at a number of different volumes. For all the data I used the lattice spacing 
of a= 0.091 fm |195j . The dependence of the nucleon mass on the length of 
the lattice looks "mild" from figure El Unfortunately, the statistical errors 
on the lightest data are too large to make any statement about finite size 
effects. 

For many years, the quality of a lattice QCD calculations was judged by 
the final value of ratio of the nucleon mass to rho mass. I now discuss some 
recent quenched QCD calculations in more detail. 

There was a large scale calculation of the quenched QCD spectrum by 
the MILC collaboration |227j that used staggered fermions. The study in- 
cluded four different lattice spacings: from 0.63 to 2.38 GeV and also 
investigated finite size effects. 

The MILC collaboration [22 7| tested 12 different chiral extrapolation 
models based on generic model for the dependence of a hadron mass Mh 
on the quark mass ntq. 

Mh = M + amy^ + brriq + cmf^ + dm^^ + emj logm^ (126) 
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Figure 16: Dependence of the mass of the nucleon on the box size in quenched 
QCD for Wilson fermions at (3 = 6.0. The data was taken from the com- 
pendium of World data in |199j . 

Only three of the fit parameters (a,b,c,d,e) were varied in a fit model. The fit 
model in equation 11261 was broad enough to include both full and quenched 
chiral perturbation theory. It is important to test more than the quenched 
chiral perturbation theory results, because it is not clear that the data is 
in the regime where the expressions are valid are valid. The simple linear 
fit, where only b was non-zero in equation I126| was not consistent with the 
data fIT7\ . 

The MILC collaboration found that the coefficient of the term from 
their fits was the wrong sign from the expectations from quenched chiral per- 
turbation theory, hence this term was not included in the their final fits. For 
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the p extrapolation the coefficient of the m^J'^ term was an order of magnitude 
lower than expected The analysis of the chiral extrapolations was complicated 
by the flavour symmetry breaking terms of staggered fermions |227j . This 
seems to be a perennial feature of the staggered fermion formalism |228j . 

The result from the MILC collaboration j227j for the ratio of the nucleon 
to rho mass ratio was m]\}/mp = 1.254 ± 0.018 ± 0.027, where the ffist error 
is statistical and the second error is systematic (in the summary figure IT7|) I 
have added the two errors in quadrature). 

Kim and Ohta |229j studied quenched QCD using staggered fermion with 
a smaller lattice spacig (a~^ = 3.63 ±0.06 GeV) and hghter quarks than were 
used by the MILC collaboration |227j . The spatial length was 2.59±0.05/m 
and their lightest quark mass was 4.5 MeV. They fitted the same chiral ex- 
trapolation models as MILC |227j did, but also had problems unambiguously 
detecting the predictions from quenched chiral perturbation theory. Kim and 
Ohta also studied some chiral extrapolation formulae, suggested by j230j . that 
looked for finite volume effects masquerading as quenched chiral logs. Kim 
and Ohta's result [221] was mjv/mp = 1.24 ± OM{stat) ± Q.Q2{sys) 

The CP-PACS collaboration found m^/m^ = 1.143 (33) (18). In the 
fits of the vector particle as a function of the quark mass, CP-PACS excluded 
the m^/^ term, but included the term The coefficient of the m^/^ term 
was found to be an order of magnitude less than the naive expectation. 

In figure IT71 I plot the results for the ratio of the nucleon and rho masses 
from some recent quenched lattice QCD calculations jH 1201 12271 1192j . I have 
selected calculations where an attempt was made to take the continuum limit. 
For comparison, in the strong coupling limit g ^ oo the hadron spectrum 
can be computed analytically. In this limit the ratio of the nucleon and rho 
masses is ^ = 1.7 HHI 

nip ' > 

The agreement with experiment for the mN/TUp ratio from quenched QCD 
is surprisingly good. Although agreement at the 10 % level may sound quite 
impressive, errors of this magnitude are too large for QCD matrix elements 
required in determining CKM matrix elements from experiment 0. The 
analysis of Booth et al. |152j using the non-analytical terms from quenched 
chiral perturbation theory estimated that the value of mpf/mp from quenched 
QCD could be as low as 1.0. This situation is not seen in the lattice results, 
however none of the quenched QCD calculations have detected the predicted 
dependence of the mass of the p meson on the quark mass. To con&m the 
quark mass dependence predicted by quenched chiral perturbation theory re- 
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Figure 17: The ratio of the nucleon mass to rho mass from several recent 
large scale quenched QCD calculations. ( jH IHHl 1223 11121 • )• 



quires new calculations with lighter quarks. The QCDSF collaboration 
obtained a coefficient of the term from their fits that was roughly in 
agreement with expectations from quenched chiral perturbation theory. 

A summary of some of the hadron spectrum from quenched QCD from 
the work by CP-PACS [96^ is shown in figure ^1 The ALPHA collaboration 
notice that the largest deviations from experiment of the quenched lattice 
results from CP-PACS are for resonant hadrons [182] if the lattice spacing 
is determined from the nucleon mass. The BGR collaboration also note 
a similar trend in their data |197j . This trend could be modelled using 
techniques similar to those used by Leinweber and Cohen |155j to study the 
p meson. It is difficult to see from figure ^1 that the hadrons with large 
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Figure 18: Final results for the light hadron spectrum from CP-PACS in 
quenched QCD 



widths have masses in quenched QCD that differ more from experiment than 
those with smaller widths. 

As discussed by Aoki [21] and the CP-PACS collaboration [HE] the value 
of the mass of the nucleon from CP-PACS [96^ and MILC ^\ differ at the 
2.5(7 level. The two calculations used different fermion formulations, each 
with a different set of potential theoretical problems, that should in principle 
produce the same results in the continuum limit. 

The main criticisms of the CP-PACS quenched study are on their treat- 
ment of the quenched chiral perturbation theory. Wittig reviews the lattice 
results for detecting the quenched chiral log 7J. The new calculations that use 
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light fermion actions with better chiral symmetry properties are disagreeing 
with the result from CP-PACS. These new calculations are done at a fixed 
(small volume) and large lattice spacings, so perhaps there are systematic 
errors in their results [7]. The CP-PACS and MILC collaborations used the 
mass of the p to determine the lattice spacing. It would have been interest- 
ing to see the dependence of the final results on using different quantities to 
determine the lattice spacing. 

Although it may be possible to do a more systematic study of quenched 
QCD at lighter quark masses than CP-PACS |196j . I am not sure that it 
is worth it. As discussed in section 13.31 the crucial non- analytic terms in 
quenched and unquenched chiral perturbation theory are very different, so 
calculations with light quarks in quenched QCD will have very little direct 
relevance for unquenched calculations. As the masses of the quarks in un- 
quenched calculations decrease we should start to see the effect of the decay 
of the hadrons. The effect of particle decay on the mass spectrum can not 
be studied in quenched QCD. 

I will now describe the results from recent unquenched lattice QCD cal- 
culations. As usual the systematic errors must first be discussed. There have 
been two recent studies of finite size effects in the nucleon from unquenched 
QCD, carried out by the JLQCD P2] and UKQCD HII]. The results are 
plotted in figure ^1 One slight concern with the nice study of finite size 
effects from JLQCD |232j is the large statistical errors on the two smaller 
volumes. Also there is some disagreement between the results from UKQCD 
and JLQCD on the 16^ volume. As stressed by the MILC collaboration [22] 
a careful control of statistical errors is required to see definitive evidence for 
the effect of the box size on the masses of hadrons. 

In figure [201 1 plot the dependence of some hadron masses on the lat- 
tice spacing from the calculations by the CP-PACS collaboration fMj. This 
should give some idea of the size of the lattice spacing errors in unquenched 
QCD. 

As discussed in section [HI the main "successes" of unquenching have oc- 
curred in the meson sector. A detailed comparison of the baryon spectrum 
with experiment is obscured by lattice spacing and finite size effects |134| 
I232j . In a preliminary analysis the MILC collaboration |233j found that 
ttim/'^p from two flavour unquenched QCD in the continuum limit was higher 
than the value in the quenched QCD, so the quenched value agreed better 
with experiment than the unquenched result. The MILC collaboration [211 
are now studying this issue using a better version of staggered fermions with 
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Figure 19: Dependence of the mass of the nucleon on the box size from 
UKQCD and JLQCD |232j . These are two flavour unquenched calcu- 
lations. 



2+1 unquenched flavours. 

In table |H1 I compare the results from lattice QCD with the old results 
from the Isgur-Karl quark model |2H4j . Although the Isgur-Karl model agrees 
better with experiment than the two lattice QCD calculations, because the 
lattice QCD calculations are based on the QCD Lagrangian, the hadron 
masses can be used to extract quark masses. This is not possible from quark 
models where there is no way of relating the constituent quark masses to the 
masses in the Lagrangian. The more interesting tests of the quark model 
occur for excited baryons. I discuss the lattice results for these hadrons in 
section 17.11 
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Figure 20: Dependence of hadron masses on the lattice spacing, from CP- 
PACS jnSl in 2 flavour unquenched QCD. Data from the finest lattice spacing 
was not used in the continuum extrapolation. 

In fact the agreement between the hadron masses from the quark model 
and experiment is actually too good. The quark model calculation of Isgur 
and Karl |2l-{4j does not include the dynamics of hadron decay. For example 
the A baryons have decay widths of around 120 MeV [Hj. It might have been 
expected that the errors in the quark model predictions for the hadron masses 
would be of the order of the decay width. Isgur and Geiger ^155| 1235] have 
developed a formalism to absorb some of the effects of quark- ant iquark loops 
into the potential. This issue is also reviewed by Capstick and Roberts |286j .. 
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Baryon 


Expt. 


Isgur-Karl MeV 


CP-PACS (quenched) 


CP-PACS (unquenched) 




940 


940 


878(25) 


1034(36) 


A 


1116 


1110 


1019(20) 


1160(32) 


S 


1193 


1190 


1117(19) 


1202(30) 




1315 


1325 


1201(17) 


1302(28) 


A 


1232 


1240 


1257(35) 


1392(58) 




1384 


1390 


1359(29) 


1488(49) 




1532 


1530 


1459(26) 


1583(44) 


Q 


1673 


1675 


1561(24) 


1680(41) 



Table 8: Comparison of the hadron masses from unquenched and quenched 
lattice QCD calculations with the resuhs from the Isgur-Karl quark 
model |234j . The unquenched data is taken from CP-PACS |134j . The 
quenched data is also from CP-PACS |196j . For the lattice results the strange 
quark mass is set by the kaon mass. 



7.1 EXCITED BARYON STATES 

Recently there has been a lot of work on the spectroscopy of excited nucleon 
states from lattice QCD. This research is mostly motivated by the experi- 
mental program at the Jefferson lab [131 CH CHI The accurate spectroscopy 
of the A^* states will be an accurate test of our understanding of the forces 
and effective degrees of freedom in hadrons |237j . Realistically, lattice QCD 
calculations may only be able to obtain one or two excited baryon states 
from a specific channel. However this is enough to (potentially) solve some 
very interesting and long standing puzzles. 

At the moment the identification of the Roper resonance from lattice QCD 
is controversial, so I will introduce some notation to prevent ambiguity. I call 
the nucleon the A^ state, the first excited nucleon with positive parity the 
A^' state, and the first excited nucleon with negative parity the A^* state. In 
the particle data table, the A^' state would be the A^(1440) and the A^* state 
would be the A^(1535). 

There is a potential artifact associated with quenched QCD for baryon 
states |238t 123^ . The mass of the rj'N state is experimentally close to the 
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Reference 


Comments 


Mn*/Mn 


Blum et al. |212! 


Domain wall 


1.49(9) 


Gockeler et al. 


Clover 


1.50(3) 


Broemmel et al. |243j 


Overlap 


1.77(7) 


Dong et al. 1210! 


Overlap 


1.67(12) 


Nemoto et al. fl^ 


anisotropic clover 


1.463(51) 



Table 9: Ratio of the mass of the parity partner of the nucleon to the nucleon 
mass from the quenched QCD. 

mass of the A^(1440) and A^(1535) states. In quenched QCD the rj' is treated 
incorrectly, so the intermediate 77' state is also incorrect. This is the ana- 
logue of the artifact in the scalar correlator found by Bardeen at al. [HS] 
discussed in section 13.11 Dong et al. |24Uj claim to have seen the correlator, 
that would be positive definite in unquenched QCD, for the N* state go neg- 
ative for pion masses below 248 MeV. If this artifact causes the correlator 
to go negative this may "confuse" fitting techniques such as the maximum 
entropy method that relies on a positive definite correlator. There is no chi- 
ral artifact in unquenched QCD. Work has started on studying the Roper 
resonance using unquenched QCD |241j . 

In table 01 1 have collected some results for the ratio of the N* mass 
to the nucleon mass from some recent quenched lattice QCD calculations. 
Only the calculation by Gockeler et al. |239j took the continuum limit. The 
experimental number corresponds to the mass of the iV(1535) divided by the 
nucleon mass about: 1.63. It is not clear what effect the 150 MeV width of 
the A^(1535) will be on the lattice result. 

The nature of the Roper resonance (iV(1440)) is still a mystery. There is 
an experimental signal for this state, but it is not clear what the quark and 
glue composition of this hadron is. On the lattice three quark interpolating 
operators are used to study the Roper state. If the mass of the Roper is not 
reproduced, then this would be evidence that additional dynamics, beyond 
three valence quarks, is important for this state. 

The quark model has problems reproducing the experimental mass of the 
Roper resonance. Using a simple harmonic oscillator potential to study the 
hadron spectrum, Isgur and Karl |242[ I245j used a oscillator quantum of 
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250 MeV. In quark model language, the N* state would have one quantum 
above the ground state, and the N' state would have two quantums above the 
ground state. This predicted ordering is opposite to the experimental masses 
of the A^(1440) and A^(1535). Capstick and Roberts review the nature of 
the Roper resonance in the context of potential models |236j . Isgur also 
discusses the problems of the Roper resonance in the quark model. |237j . 
The predictions of the quark model for the lowest excitations of the nucleon 
improve if a more realistic potential is used, and the mixing between states 
is taken into account |236| 1237] . 

There are predictions from flux tube and bag models that the lightest 
hybrid baryon (three quarks with excited glue) |24ti| I247j is = l"*" with 
a mass in the region 1.5 to 1.9 GeV, hence close to the mass of the Roper 
resonance. 

Sasaki, Blum, and Ohta |242j studied the first excited state of the nucleon 
at a fixed lattice spacing of 0.1 fm and with a physical length of the lattice as 
1.7 fm. Their calculations were done in quenched QCD with quark mass in 
the range Mps/My 0.59 to 0.9. The excited state masses were extracted using 
a variational technique with two basis states that were different interpolating 
operators for the nucleon. Sasaki et al |242j could only obtain a signal for the 
N' state for pion masses above 600 MeV. The mass of the A^' state was larger 
than the A^* state. If the variational technique was not biased by truncation 
of the sum of excited states, then the calculation of Sasaki |242j should be 
able to resolve the ground and first excited states. For the negative parity 
states, the masses obtained from the ground and first excited state were 
degenerate. This could be interpreted as the variational technique not being 
able to resolve two states. Experimentally the lightest A^* states are the 
N{1535) and Ar(1650). 

Melnitchouk et al. |248j have also studied the spectrum the masses of the 
N, N*, and N' states. Their raw data for the splitting between the N and 
N* states is consistent with other groups |242t I239j . The mass for the A^' 
state was much higher than the mass of the experimental Roper resonance. 

Dong at al. |240j claim agreement between the mass of the A^' state from 
their calculation and the mass of the Roper resonance from experiment. The 
calculation was done with the lattice spacing of 0.2 fm and a physical box 
lengths of 2.4 and 3.2 fm. The calculations used very light valence quarks (for 
the lattice anyway). The lightest pion mass was 180 MeV. The excited state 
masses were extracted using the constrained curve fitting method developed 
from the proposal by Lepage at al. [Ul]. Dong at al. |240j reported that the 
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mass of the N' state started to decrease rapidly with pion masses below 400 
MeV. The mass of the A^' state was less than that of the N* state for pion 
mass of about 220 MeV. The mass quoted for the A^' state is 1462(157) MeV. 
The large error bars on this result means that a 2 a statistical fluctuation 
would give 1776 MeV. 

Sasaki |249j has reported a study of the finite size effects on the mass of 
the excited state of the nucleon. The excited state of the nucleon was studied 
at the fixed lattice spacing of 0.09 fm. Three physical box sizes were used: 
1.5 fm, 2.2 fm and 3.0 fm. The mass of the N' state was extracted using the 
maximum entropy method. The A^' and N states had finite size effects that 
increased as the light quark masses were reduced. The final result for the N' 
state looks as though it would extrapolate to the experimental value from 
pion masses around 600 MeV. No dramatic decrease in the mass of the N' 
is required. The picture from the calculation of Sasaki |249j disagrees with 
that from Dong et al. |240j . 

Broemmel et al |243j have tried to study the Roper resonance using lattice 
QCD. They used an overlap-Dirac fermion operator with a lattice spacing of 
0.15 fm and two physical lattice sizes of 1.8 and 2.4 fm. To study excited 
states they used a variational technique based on three interpolating opera- 
tors. They could get a signal for the nucleon with pion masses as low as 220 
MeV. Unfortunately, they could extract a signal for excited nucleon states 
with pion masses above about 550 MeV. This stopped them being able to 
confirm the mass dependence of the A^' state claimed by Dong et al. |240j . 
This analysis did not include the effect of the quenched chiral artifact in this 
channel. Broemmel et al |243] extracted both the ground and first excited 
state of the negative parity nucleon channel. 

From the very interesting recent studies |242t I240( I249( I243j of the Roper 
resonance on the lattice it is clear that careful attention will have to be paid 
to the systematic errors. Pion masses below 200 MeV may be required. It 
would be good to have variational calculations with wider sets of interpolating 
operators as a check on the various Bayesian based fitting techniques. For 
example it would be reasonable for "fuzzed" nucleon operators jlH] to couple 
more to hybrid baryon states because they contain more glue. It is a high 
priority that other groups try to reproduce the results of Dong at al. |240j . 
and Sasaki [249J. 

One particular concern with getting the mass of the A^(1440) is that its 
width is 380 MeV. The decay of A^(1440) will effect its mass. This is also 
a difficulty for potential model calculations. Capstick and Roberts |236j 
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estimate that ignoring the width of the A^(1440) causes an uncertainty of the 
order of 100 MeV on the mass within the potential model framework. 

Although I have focused on the excited states of the nucleon, there is 
growing body of work on the parity partners of other baryons |244[ I248j . 

8 ELECTROMAGNETIC EFFECTS 

The majority of lattice QCD calculations do not incorporate the effect of 
the electromagnetic fields in hadron mass calculations. This is reasonable, 
because the dominant interaction for hadron masses is the strong force. I 
review the work done on including electromagnetic fields in lattice QCD 
calculations. 

The theory of QED has been studied by many groups using lattice tech- 
niques. The formalism is similar to that discussed in section |21 except that 
the gauge group is U(l). However, there are some conceptual differences, 
because the U(l) gauge theory is not asymptotically free. Issues relating 
to the non-asymptotically free nature of QED have been studied in a non- 
perturbative way on the lattice |2501 1251j . Also there are compact and non- 
compact versions of the lattice U{1) theory. In non-compact QED the gauge 
fields A^^ take to the range — oo to oo. 

To study electromagnetic effects on the hadron spectrum the electromag- 
netic fields have been quenched. The dynamics of the sea quarks have not 
been included in the generation of the U(l) gauge fields. A comparison be- 
tween the use of background fields in sum rules and lattice QCD has been 
made by Burkardt et al. |252j . 

The most "comprehensive" study of the effect of electromagnetism on 
the masses of hadrons has been performed by Duncan et al. |253t 1254] . They 
used a non-compact version of QED The gauge fields were generated using 
the action 

Sem = ^Yl EiD.A^i^) - D^A.ix)) (127) 

where e is the electromagnetic charge. The fields were subject to 

the linear Coulomb condition. The fields were promoted to compact fields 
U{x)fj_ = e='=*''"^''(^). This field coupled to a quark field with charge ±ge. 

Currently, lattice QCD calculations are usually not accurate to 10 MeV, 
the order of magnitude of electromagnetic effects on the masses of light 
hadrons. To increase the mass splitting Duncan et al. |25H| I254j used large 
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Mass splitting 


Raw Lattice QCD (MeV) 


Corrected MeV 


Experiment (MeV) 


m7r+ — Tn^_ 


4.9(3) 


5.2(3) 


4.594) 


rrin - rrip 


2.83(56) 


1.55(56) 


1.293 


m^o — ms+ 


3.43(39) 


2.47(39) 


3.27 


ms- - m^o 


4.04(36) 


4.63(36) 


4.81 


m^- — m=o 


4.72(24) 


5.68(24) 


6.48 



Table 10: Electromagnetic mass splittings from Duncan at al. |253| I254j 

charges (2 to 6 times the physical values) and then matched onto chiral per- 
turbation theory that included the photon field. 

The chiral extrapolation fit model used for pseudoscalars with electro- 
magnetic fields 

mp = A{eq, e-q) + mgB{eq, e-g) + mqB{eq, e-q) (128) 

where Cg, e-q {rUq, m-q) are the (masses) of the quark and anti-quark 

The calculations of Duncan et al. |253| I254j were done with one coarse 
lattice spacing ~ 1.15 GeV with a box size of 2 fm. Some of the results 
for the mass splittings from Duncan et al. are in table ^1 I have included 
both the data and the corrected results, where theoretical expressions were 
used to correct for finite volume effects. 

There has some work on computing electromagnetic polarizabilities |255| 
I256t 1257] from lattice QCD. The electric {E) and magnetic {B) polarizabili- 
ties measure the interaction of a hadron with constant electromagnetic fields. 
Under the electromagnetic field the mass of the hadron is shifted by 5m. 

Sm = ~aE^ + (129) 

The a and (3 quantities should be computable from QCD. Holstein [258] re- 
views the theory and experiments behind the nucleon polarizabilities. The 
experimental values for a and (3 are extracted from Compton scattering ex- 
periments (see [259J for example3). A comparison of the results from lattice 
QCD to models and experiment can be found in two recent papers |256[I257] . 
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In contrast to work of Duncan et al. the formalism used for electromag- 
netic polarizations |255j uses static electromagnetic fields. The SU (3) gauge 
fields are modified by multiplying 

Uiix) ^ e""'^'''Ui{x) (130) 

where X4 is the Euclidean time variable and E is the constant electric field. 
The phase factor in equation 11301 can be linearised. Smit and Vink have 
described how to put a constant magnetic field on the lattice |260j . 

There are speculations that in very strong magnetic fields [B > 5 x 10^^ 
T), the proton will become unstable to the decay to neutrons |261j . Magnetic 
fields of this intensity may be realised in the universe j262j . The original 
estimate |261j of the instability of the proton was done in the quark model. 
In an attempt to remove some of the uncertainty in the hadronic calculation, 
Rubinstein et al. |253j used lattice QCD to study the dependence of the 
masses of the proton and neutron on the magnetic field. 

Some early lattice calculations included the magnetic fields to look at 
the magnetic moments of hadrons j2641 1265j . However, it is best to calculate 
magnetic moments from form factors |266[ 12671 1268j . so electromagnetic fields 
are no longer used. 

Electric fields were used in the first (unsuccessfully) attempts to compute 
the electric dipole moment of the neutron |269t 1270] . The calculation of the 
neutron electric dipole moment has recently been reformulated [271j in a way 
that does not require the use of electromagnetic fields. 

9 INSIGHT FROM LATTICE QCD CAL- 
CULATIONS 

The start of the book Numerical Methods for Scientists and Engineers by 
Hamming has the immortal phrase: "the purpose of computing is insight 
not numbers." In the previous sections I have described how lattice QCD 
is used to compute the masses of hadrons. This may give the impression 
that lattice QCD is essentially just a black box that produces the masses of 
hadrons without any insight into the physical mechanisms or relevant degrees 
of freedom. In this section I hope to show that lattice QCD can also help 
to explain the physical mechanisms behind the hadron mass spectrum. A 
very good overview of the type of insight wanted from hadronic physics, that 
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contrasts the hadron spectroscopy approach to the study of confinement with 
the results from DIS type studies, is the paper by Capstick et ah |272j . Note 
however for the B physics experimental program, high precision numbers 
with reliable error bars are required to look for evidence for physics beyond 
the standard model of particle physics |273j . 

Isgur's motivation |237j for studying the N* particles is based on trying to 
understand the important degrees of freedom that describe low energy QCD. 
Increasingly, lattice QCD calculations are being used to provide insight into 
the dynamics of QCD. Some of Isgur's last papers |274| I275t llU5t llU6j were 
devoted to using lattice QCD to validate the quark model picture of hadronic 
physics. 

At first sight the model building approach to studying QCD appears to 
give more insight into the dynamics of QCD. It is usually quite easy to study 
the effect of adding new interactions to the model. There are many mod- 
els for QCD interactions: quark model j236j . instanton liquid model |276j . 
or bag model |277j . The main problem with the model based approach to 
hadron spectroscopy is that is very difficult to judge whether the assumptions 
in the models are valid. The models are only believed when they provide a 
reasonable description of "most" experimental data with a good x^, but this 
does not prove that they are correct. Different models of the QCD dynamics 
can be based on very different physical pictures, but may give equally valid 
descriptions of experimental data. For example the physical assumptions 
behind the bag model seem to be very different to the assumptions behind 
the instanton liquid model |276j . The major advantage of lattice QCD cal- 
culations is that the theory can be mutilated in a controllable way, so the 
physical mechanism underlying a process can be studied. 

For example, the question about what mechanism in QCD causes a lin- 
early rising potential for heavy quarks at intermediate distances is not some- 
thing that can be answered by the quark model. Greensite |i278j reviews the 
work on studying "confinement" using lattice QCD). 

In this section I will focus on various attempts to explain the value of 
the mass splitting between the nucleon and delta. According to the quark 
model [U], the masses of the nucleon and delta are split by a spin-spin hyper- 
fine term: 

Hhf = -asMj2{aXaUaXa)j (131) 

i>j 

where the sum runs over the constituent quarks and is the set of SU(3) 
unitary spin matrices (a runs from 1 to 8). In perturbative QCD a term of the 
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Figure 21: Edinburgh plot from CP-PACS from quenched QCD. The 
lattice data is compared against the quark model of Ono |281j . 



form 113 ll is naturally generated, but it is not clear whether the equation 11311 
is relevant for light relativistic quarks. 

Lattice QCD calculations can produce more masses than experiment, so 
the hyperfine interaction can be tested |279t l28Uj . A pictorial measure of the 
quark mass dependence of the masses of the octet and decuplet is given by 
the "Edinburgh" plot from CP-PACS [Hn] is in figure |^ The continuous 
curve is from a quark model by Ono |281j that uses: 

S-S- 



Mmeson = + TU, + ^^T. (132) 

The agreement between the lattice data and Ono's model is reasonable in fig- 
ure [1221 The lattice QCD data from CP-PACS ^ in figure Ell is almost pre- 
cise enough to show deviations from the model by Ono |281j (equation I132|l . 
In the instanton liquid model |l()lj . the vacuum is made up of a liquid 
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of interacting instantons. There has been a lot of work on comparing the 
instanton hquid model against lattice QCD. The basic idea is to cool the 
gauge configuration [21]. This removes the perturbative part of the gauge 
field and leaves the classical configurations, that can be compared against 
the predictions of the instanton liquid model. The cooling procedure is es- 
sentially a way of smoothing the perturbative noise from the configurations. 
This perturbative noise presumably has something to do with the one gluon 
exchange term. 

In the first study of this on the lattice |282j . the masses of the nucleon, p, 
A and vr were measured in the usual way. The gauge configurations were then 
cooled and the simple hadron spectrum was measured again. The masses for 
the p IT and nucleon particles were qualitatively the same before and after 
cooling. As the cooling does not effect the instantons this suggested that the 
mass splittings are largely due to instantons. The mass splitting between 
the nucleon and A was reduced by smoothing. Chu et al. j282j claimed 
that this was due to the problems with extracting the mass of the A from 
the lattice data. Later work on this issue |283| llUUj has not returned to the 
effect of instantons on the mass splitting between the nucleon and A. Rosner 
and Shuryak [2S^ have shown that simple instanton interactions can give a 
reasonable representation of some baryon mass splittings. 

In a series of papers, Glozman |285j and collaborators have argued that a 
interaction based on Goldstone boson exchange between constituent quarks 
gives a better description of the mass spectroscopy of baryons, than interac- 
tions of the form in equation 11311 

This is not the place for a detailed review of the case for and against a 
interaction based on the exchange of Goldstone Bosons (GBE). For a critique 
of GBE you can look at the paper by Isgur 286 and the review article by 
Capstick and Roberts [236^. Here I will just discuss the evidence for GBE 
from lattice QCD. 

The Kentucky group |287j have introduced valence QCD, a mutilated 
version of lattice QCD, that omits "Z" graphs from the formalism. The aim 
was to study the foundation of the quark model. In quenched QCD the 
sea quark loops are omitted, however there are still higher Fock states from 
intermediate "Z" states. Figure |221 shows that a quark going backwards in 
time can be interpreted as meson state. 

The lattice version of valence QCD is the Wilson fermion action (equa- 
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Figure 22: Z graph for a quark. 
tionlH)) with the backwards hopping terms removed. 

X i=l 

+ ^x^x - i^i^x+iili + l)Ul{x)i)^ + - 7t)^x) 

The Kentucky group |287j studied valence QCD in a lattice calculation 
aX (3 = 6.0 with a volume of 16^ 24, and sample size of 100. The Kentucky 
group |287j also used valence QCD to study form factors and matrix elements, 
but I just will focus on their results for the hadron spectrum. 

In figure 1221 I show a comparison of some hadron masses (from the com- 
ment by Isgur j274j ) The main conclusion from figure ESI is that the hyperfine 
splittings seem to have been reduced for light hadrons. 

Valence QCD still has the physics of gluon exchange, so the near degen- 
eracy of the nucleon and A, suggests that the "Z" graphs are the important 
mechanism behind the nucleon-delta mass splitting. For heavy hadrons the 
results from VQCD are not suppressed relative to the those from quenched 
QCD. For example, the VQCD prediction for the mass splitting between 
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Figure 23: Spectrum of valence (vQCD) and quenched (qQCD) QCD 
from ^I^. 

and B is 45.8 ± 0.4 MeV, compared to the experimental value of 46 MeV. So 
the one gluon exchange is important for heavy quarks (beyond charm), but 
less important for hght quarks. The VQCD calculations have much reduced 
hadron masses. For example, the nucleon mass in VQCD is reduced by 700 
MeV from quenched QCD. The reduction of the vector meson mass was 537 
MeV. The Kentucky group |287j attributed the reduction in hadron masses 
to reduction of the dynamical or constituent quark mass due to the omission 
of Z graphs. 

Isgur j274j criticised the conclusions from the study of lattice VQCD. In 
particular, Isgur notes that the one boson exchange interactions, based on 
Z graphs, operates between two quarks, so should not effect mesons made 
from a quark and anti-quark. So the near degeneracy between of the p and vr 
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masses in figure ESI would not be expected from tlie suppression of Z graphs. 
Also Isgur |274j noted that the hadron spectrum from valence QCD was 
radically different from experiment (and quenched QCD). This might cause 
additional problems, if for example the wave-functions were very different for 
valence QCD compared to the real world, then this would effect the hyperfine 
splittings. The Kentucky group did study matrix elements that are related 
to the wave function of the states, but did not see any problems |288j . 

The current lattice calculations are currently not able to determine the 
correct mechanism for the mass splitting between the A and nucleon. Much of 
the lattice work has concentrated on validating one particular model, rather 
than also falsifying competing models. This type of physics is necessarily 
qualitative. For example, in the spirit of theory mutilation, ideally only one 
piece of physics must be removed at one time. It is not clear what effect that 
valence QCD has on the instanton structure of the vacuum. The modification 
of the Wilson fermion operator in equation 11331 will also effect the zero mode 
structure that is of crucial importance to instanton inspired models. The 
operator in equation 11331 does not obey 

Mt = 75M75 (133) 

for the Wilson operator M. There is a generalisation of equation 11331 in va- 
lence QCD. There has been some work on trying to "disprove" the instanton 
liquid model by studying lattice QCD gauge configurations flU5]. Edwards 
reviews the work by many groups on this |289j . 

Lattice QCD simulations have been used to test other assumptions made 
in models of the QCD dynamics. For example, there are some models of 
hadronic structure that are based on diquarks |29(Jj . A critical assumption 
in diquark models is that two quarks actually do cluster to form a diquark. 
This assumption has been tested in a lattice gauge theory calculations by 
the Bielefeld group |291j . where they found no deeply bound diquark state 
in Landau gauge. Leinweber |292j has claimed that lattice QCD data on the 
charge radii of hadrons provides evidence against scalar diquark clustering. 
As a test of the MIT bag model |293j and Skyrme model |294j . density 
operators from the models were compared against results from quenched 
QCD. 

The large A*";, limit of QCD provides much insight into QCD (see |295j for 
a review). Teper |296j and collaborators have studied the glueball masses in 
the large Nc limit. Lattice QCD calculations can be done with any gauge 
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group. Teper |296j et al. studied the glueball masses for Nc = 2,3,4,5. This 
allowed them to estimate the size of the corrections to the N^. ^ oo limit. 
As the 0{1/N^) corrections are small, it was important to control both the 
finite volume and lattice spacing errors. 



10 WHAT LATTICE QCD IS NOT GOOD 
AT 

In the previous sections I have implicitly assumed that all the hadrons are 
stable. In the real world, most hadrons are unstable to strong decays. For 
example the p meson has a mass of 770 MeV and a decay width of 150 
MeV. Most lattice practitioners never "worry" about the p's decay width 
(in public at least). The determination of a hadron's mass from experi- 
ment is inextricably linked to the determination of the decay width. This is 
perhaps most dramatically demonstrated by the problem of missing baryon 
resonances [297J. The quark model predicted more excited baryons j245] 
than were actually seen in A^vr reactions. It was claimed that the additional 
states were not seen because they coupled very weakly to the Ntt channel. 
The quark model did predict that the missing baryons might be seen in Ninr 
reactions |297j . There are experiments at the Jefferson lab that are trying to 
detect these "missing resonances" [H] . 

Lattice QCD calculations have to be done in Euclidean space for conver- 
gence of the path integral |35j. This implies that the amplitudes and masses 
from lattice calculations are real. This makes the study of resonances non- 
trivial, because decay widths are inherently complex quantities. This is also 
a problem for calculations with a finite chemical potential, although there 
has been some progress in this area |298j . 

I do not discuss the very elegant formalism of Liischer |299j for studying 
unstable particles in a finite volume. The formalism and results from the 
scattering formalism are reviewed by Fiebig and Markum in this volume |300j . 
In this section, I would like to describe the possible implications for mass 
determinations from standard correlation functions. 

The momentum is quantised on a lattice of length L and periodic bound- 
ary conditions. The momentum of mesons can only take values: 

2?™ , 

Pn = —r 134 

aL 
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where n is an integer between and L — 1. For a typical lattice, L = 16 and 
= 2.0 GeV, so the quantum of momentum is 0.79 GeV. The quantisation 
of momentum makes the coupling of the state to the scattering states different 
to that in the continuum. This "feature" has been used to advantage by 
Liischer |299j in his formalism. 

The quantisation of momentum has important consequences for mesons 
that decay via P wave decays such as the p meson. A p meson at rest can 
only decay to two pions with momentum p and —p. The decay threshold 
is 2^m1 + (^)^. The quantisation of momentum on the lattice does not 
effect the threshold for S"— wave decays, An example of a S-wave decay is the 
decay of the O^"*" meson to pairs of mesons. The quarks in current lattice 
calculations are almost light enough for the strong decay of the flavour singlet 
0++ fSUT]. 

Hadron masses are extracted from lattice QCD calculations using two 
point correlators f equation I23|) . However the use of equation mav not 
be appropriate for hadrons that can decay. The most naive modification of 
the lattice QCD formalism caused by the introduction of decay widths is the 
replacement 

m — > m + — (135) 

where F is the decay width. This modifies equation 1221 

c(t) = aoe-™»*e-^ + ■ ■ ■ (136) 

that is used to extract the masses from correlators. There are a number of 
problems with equation 1136^ so a more thorough derivation is required. 

I review the work by Michael |302j (see also the text book by Brown |303j 
) on the effect of decays on two point functions. I will consider two scalar 
fields (a, vr) interacting with the interaction crTrvr The mass of the vr particle 
is /i. To study the implications of particle decay on the two point correlator, 
consider the renormalisation of the propagator of the o particle in Euclidean 
space 

PB\p)=P^ + m^ (137) 

The effect of the interaction of the tt particle with the a particle renormalises 
the a propagator. 

p-\p)=p'' + m^-X{p'') (138) 
where X is the self energy. 
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Masses are extracted from lattice QCD calculations using the time sliced 
propagator: 

G{t) = — / cipoe^^"*P(po,Q) (139) 
= - / dEe-^'p{E) (140) 



where 

If the pole around E ~ —M is neglected then the expression for the 
spectral density is 

In the limit 7 < (M - 2/i), (M - 2/i)t > 1 

G(t) = — e-*" cos(7t) + J- (143) 

In figure |^ I plot separately the log of the first and second terms in 
equation 11431 using the parameters: (M = 0.5,7 = 0.05, /i = 0.1). The 
breakdown of a pure exponential decay can be seen for the second term. 

There is no evidence for the breakdown of the simple exponential fit 
model in current lattice calculations (apart from the effects of chiral artifacts 
in quenched QCD |HS1)- This effect may become apparent as the sea quark 
masses are reduced to where particle decay is energetically allowed. 

The maximum entropy approach |304t 17^ to extracting masses from cor- 
relators produces an estimate of the spectral density, so in principle could be 
used to extract the decay widths of particles. It is not clear to me whether a 
decay width obtained from an analysis based on maximum entropy method 
would be physical. Yamazaki and Ishizuka |305j have recently compared the 
maximum entropy approach to studying unstable particles to the method ad- 
vocated by Liischer 1299) in a model theory. Yamazaki and Ishizuka claimed 
good agreement between the two methods. 

Another way of looking for "evidence" of resonant behaviour is to look 
for peculiarities in the quark mass dependence of the hadron masses. One of 
the first applications of this idea was by DeGrand who studied the effect of 
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Figure 24: Resonant correlators from equation 11431 

p decay on the quark mass dependence of the rho mass jl54j . The mixing of 
the p correlator with tttt states at the nonzero momentum makes the mass 
dependence of the p correlator more complex than described by the theory in 
section Eini The MILC collaboration jSH] have recently claimed to see some 
evidence for the decay of the qq (non-singlet 0^^) meson by comparing the 
mass dependence of the oq particle in unquenched and quenched QCD. 

The MILC collaboration PHj tried to look for the evidence of p meson 
decay by studying the dispersion relation of the rho particle. The momentum 
can be injected into particle correlators, so that the hadron masses can be 
computed at nonzero momentum. The signal to noise gets worse, as the 
momentum increases, so typically the dispersion relation is only known for 
a few values of momentum. On the lattice the dispersion relation is of the 
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form |3(J(jj 

sinh^ E = sinh^ m + sin^ p (144) 

The "improvement" program described in the appendix, aims to make the 
lattice dispersion relation closer to the continuum one. Some clever lattice 
theorists call this computing "the speed of light" |lU8j . 

The relationship in equation 11441 predicts that the dependence of the 
hadron mass as a function of momentum. Consider an interpolating operator 
for the p meson: ip'-fiip. If a unit of momentum is injected in the z direction, 
then the iplz'^P operator will couple to the p with one unit of moment, as 
well as two pions with unit momentum. The two states will mix and the 
masses will be modified. However, the operators: 'ip'jyil', ip-^x^J do not couple 
to the two pion states, essentially because pyrvr interaction is zero for these 
kinematics. A signal for rho decay would be a different mass from the p's 
polarised perpendicular and parallel to the momentum of the state. The 
MILC collaboration did not see this effect, perhaps because of artifacts with 
the type of fermions used. The UKQCD collaboration have recently claimed 
to see evidence for p decay via this mechanism |3U7j that was consistent with 
the results from other methods. 

In this review I have focused on just computing the masses of hadrons. 
The computation of decay widths from first principles is the holy grail of 
hadron spectroscopy. There have been very few attempts at the calculation 
of the decay widths from lattice QCD. The hadronic coupling constants are 
computed from matrix elements, rather than trying to fit expressions like 
equation 11431 to data. Additional correlators over the standard two point 
correlators are required to be calculated. There have been some attempts to 
compute the rho to two pion coupling constant |3U8| l3U9t I31U1 l5U7j . Perhaps 
the most famous calculation decay widths is the calculation of partial widths 
for the decay of 0+^ glueball were calculated by the GFll group |311j . A 
clarification of the formalism for studying decay widths and mixing in lattice 
QCD calculations is described in |3U21 lUU] . 

11 CONCLUSIONS 

The computation of the light hadron mass spectrum is essential to checking 
the validity of lattice QCD techniques. For example, to reliably extract quark 
masses requires a precise and consistent calculation of the masses of the light 
mesons. A precise determination of the light hadron mass spectrum would 
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unequivocally demonstrate that non-perturbative quantities can be extracted 
from a physically relevant quantum field theory. 

The recent large scale quenched lattice QCD calculations of the spectrum 
of the lightest hadrons from CP-PACS jHEl I196j are the benchmark for future 
calculations. After nearly twenty years of work the deviations of the predic- 
tions of quenched QCD from experiment have been quantified. The next 
generation of large scale quenched QCD calculations (if they are worth doing 
at all) will probably use the new fermion operators, such as overlap-Dirac 
operator, to push to lighter quark masses below Mps/My of 0.4. The light 
quark mass region of quenched QCD is full of pathologies, such as quenched 
chiral logs, that will be fun to study theoretically, but of limited or no rele- 
vance to experiment. Quenched QCD calculations will still be of some value 
for many interesting quantities where the uncertainty is larger than the in- 
herent error of quenched QCD. 

The main challenges in unquenched hadron spectrum calculations is de- 
termining the lattice spacing dependence and reducing the size of the quark 
masses used in the calculations. The unquenched calculations of the CP- 
PACS collaboration |134j have shown that at least three lattice spacings will 
be needed to obtain high "quality" results. Lattice QCD calculations with 
2+1 fiavours of quarks are starting to produce important results The 
effective field theory community is increasingly doing calculations specifically 
to analyse data from lattice QCD calculations. 

Motivated by the "new nuclear physics" experimental programs at facil- 
ities, such as Jefferson lab dH [T3], lattice calculations are starting to 
be used to study interesting particles such as the N*s. The use of more so- 
phisticated interpolating operators and more advanced statistical techniques 
(such as maximum entropy techniques) may allow some information to be 
obtained on some of the lowest excited states of hadrons. The determination 
of the masses of excited states from lattice QCD would be a big step forward 
for hadron spectroscopy, if this was indeed possible. 

The theory of hadronic physics is the ultimate postmodern playground, 
as it sometimes seems that the use of a particular model for hadronic physics, 
from a mutual incompatible set of possibilities, is almost a matter of personal 
preference |3 12j . One advantage of lattice QCD calculations is that they pro- 
vide qualitative information about physical mechanisms, that is not directly 
accessible from experiment. If as Shuryak |276j claims, the physical picture 
behind the bag model and the instanton liquid model are different, then only 
one picture is correct, so one of them must be discarded. I hope that qualita- 



88 



tive lattice QCD calculations can help simplify the theories behind hadronic 
spectroscopy by ruling out the underlying pictures behind certain classes of 
models. The aim of simplifying the theory of hadronic physics is a stated 
aim of the current experimental program |2371 1272j . 

There are number of interesting "challenges" for hadron spectrum calcu- 
lations beyond the critical task of reducing the errors in lattice QCD calcu- 
lations: 

• Can the resonant nature of the p meson be determined from lattice 
QCD? 

• Can improved lattice calculations determine the structure of the Roper 
resonance (first excited state of the nucleon). For over thirty five years, 
there have been many speculations on the nature of the Roper reso- 
nance. Can lattice QCD close this issue? 

• Can the physical mechanism behind the hyperfine splittings in mesons 
and baryons be determined from lattice QCD? 

• Can lattice QCD calculations be used to simplify hadronic physics by 
ruling out (or perhaps even validating) the bag model? 
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A TECHNICAL DETAILS 

To outsiders (physicists who live in continuum), the lattice gauge theory 
community most seem like a very inward looking bunch. A large fraction of 
research in lattice gauge is on improving the methods used in lattice calcu- 
lations. Improving the methods used in lattice calculations leads to smaller 
and more believable error bars and hence is a good thing! Currently, the 
biggest improvements in the methodology of lattice QCD calculations are 
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coming from "better" lattice representations of the continuum Dirac oper- 
ator. The "new" lattice representations of the Dirac operator have either 
reduced lattice spacing dependence or a better chiral symmetry. 

The importance of the dependence of the results on the lattice spacing has 
been stressed through out this review. As the results of lattice calculations 
are extrapolated to the continuum, the calculations would be more precise if 
the lattice spacing dependence of quantities was weak. The computational 
cost of reducing the lattice spacing used in lattice QCD calculations from 
equation HO is very large, hence it is advantageous to use coarser lattice 
spacings 107, 313,. 

A standard technique from numerical analysis is to use derivatives that 
are closer approximations to the continuum derivatives. For example the 
lattice derivative in equation 11461 should be more accurate with a larger 
lattice spacing than derivative in equation | 



/(x + a)-/(x-a) ^ j/^ (145) 
2a ax 

4 /(x + g) - /(x - g) /(x + 2a) - f{x - 2a) df . 

3^ 2a IQa ^ dx ^ ' ^ ' 

However in a quantum field theory there are additional complications, such 
as the operators in equation 11461 mixing under renormalisation. 

There is a formalism due to Symanzik |314[ I315j called improvement 
where new terms are added to the lattice action that cancel 0(g) terms 
(irrelevant operators) in a way that is consistent with quantum field theory. 
The required terms in the improved Lagrangian can be simplified by the use 
of field redefinitions in the path integral |316j . A very elegant numerical 
procedure to improve the Wilson fermion action has been developed by the 
ALPHA collaboration (see the review [189] by Liischer). 

The Wilson fermion operator in equation El differs from the continuum 
Lagrangian by 0(g) terms. The improvement scheme used in most lattice 
QCD calculations with Wilson fermions is called clover improvement. The 
clover term |316j is added to the Wilson fermion operator in equation |H1 

^clover ^ ^ ^(^a.^F,^^.) (147) 

where F^^ is the lattice field strength tensor. 

If the csw coefficient computed in perturbation theory is used then the 
errors in the results from the lattice calculation are 0{ag^). The ALPHA 
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collaboration |317j have computed csw to all orders in g'^ using a numerical 
technique. The result for csw from ALPHA is 



Csw — 



1 - 0.656^2 _ 0.152^^ - 0.05%6 
1 - 0.922^2 



(148) 



for < 5^ < 1, where g is the coupling. The estimate of csw, by ALPHA 
collaboration, agrees with the one loop perturbation theory for g < 1/2. 

The clover improvement program for Wilson fermions has had many prac- 
tical successes. Unfortunately, it is computationally very costly to reach 
light quark masses in quenched or unquenched lattice calculations that use 
the clover fermion operator |13Uj . Hence, attention has focused on also im- 
proving the eigenvalue spectrum of the lattice representation of the Dirac 
operator. 

The design of fermion operators on the lattice has a deep connection with 
chiral symmetry and the global chiral anomaly. The theoretical complications 
with transcribing the Dirac operator to the lattice are reviewed in many 
places pm [T7] . Our understanding of chiral symmetry on the lattice has 
recently increased by the rediscovery of the Ginsparg- Wilson relation |818j 



where M is the fermion operator in equation |H1 at zero mass. Equation 11491 
smoothly matches onto the chiral symmetry equation in the continuum as 
the lattice spacing is taken to zero. Lattice fermion operators that obey 
the Ginsparg- Wilson relation (equation I149|) have a form of lattice chiral 
symmetry |319j . Explicit solutions, such as overlap-Dirac |320j or perfect 
actions |321l|. to equation 11491 are known. Actions that obey the Ginsparg- 
Wilson relation are increasingly being used for quenched QCD calculations |322j . 
Domain Wall actions, that can loosely be thought of as being approxi- 
mate solutions to the Ginsparg- Wilson relation, have been used in calcu- 
lations |323t I324j of the matrix elements for the e'/ e. 

The main downside of fermion operators that obey the Ginsparg- Wilson 
relation is that they are computationally expensive. In a review of the lit- 
erature, Jansen argues |325j that overlap-Dirac type operators are roughly 
100 times more expensive computationally than calculations with standard 
Wilson fermions. The development of new algorithms should reduce this 
difference in computational cost |197j . The various versions of the overlap- 
Dirac operator are cheap enough to use for quenched calculations. I speculate 



M75 + 75M = aM75M 



(149) 
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that there will be an increasing trend to use overlap-Dirac type operators for 
quenched calculations. It will be some time before unquenched calculations, 
that are useful for phenomenology, are performed with overlap-Dirac opera- 
tors. Unquenched calculations of QCD with domain wall fermions have just 
started [326J. 

A more pragmatic development in the design of light fermion actions 
is the development of improved staggered fermion actions [201 1327j . This 
class of action is being used for unquenched lattice QCD calculations with 
very light quarks (see table Ej) by the MILC collaboration. The problem 
with standard Kogut-Susskind quarks was that the formalism broke flavour 
symmetry. So numerical calculations usually had fifteen pions split by a 
considerable amount from the Goldstone boson pion. The new variants of 
fermion operators in the staggered formulations have much reduced flavour 
symmetry breaking. The improved staggered quark formalism is quite ugly 
compared to actions that are solutions of the Ginsparg- Wilson relation, but 
lattice QCD is a pragmatic subject and utility wins out over beauty. It is 
not understood why calculations using improved staggered quarks are much 
faster ||82j than calculations using Wilson fermions jSlj . 
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